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ABSTRACT 


The fraction of the light in a stellar image that passes the slit of a spectrograph is 
calculated, and the limitations on any increase of this fraction by changes in the lens or 
mirror system are discussed. 

It is shown that most of this light may be made to pass the slit by means of a device, 
here described, which slices up the image into a series of strips of a width equal to that 
of the slit and which arranges them end to end along the slit. The wide spectrum thus 
produced can be narrowed down by means of a cylindrical lens, thereby increasing the 
intensity of illumination at the photographic plate by factors ranging up to 10 or more. 

The aberrations of the cylindrical lens used to narrow the spectrum are discussed. 


I. EFFICIENCY OF A STELLAR SPECTROGRAPH IN THE 
USUAL ARRANGEMENT 
The image of a star or other object, subtending an angle a, has 
a diameter 
a = aD,F, (1) 


at the focus of a telescope of aperture D, and focal ratio F,. If the 
spectrum of this stellar image is to be studied by a spectrograph 
whose collimator and camera lens have effective focal ratios F, and 
F,, the slit width must be 


in 22 (2) 
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in which A is the width of the slit image on the plate. If the full 
resolving-power of the spectrograph is to be attained, A must be 
made equal to, or less than, the resolving-power of the photographic 
plate, i.e., 20-40 wu. Assuming uniform illumination over a circu- 
lar image, the fraction (R) of the light in the image passed by the 
slit is, then, approximately 

it 4ab—4F A ( 


T a? wh aD,.F | ; 


wn 
— 


However, if all the light entering the spectrograph is to pass the 
collimator lens, it is necessary that F. < F;. The maximum ef- 
ficiency of the slit therefore is 


- 1... (4) 
waD.F, 
With the large-aperture telescopes that are now becoming avail- 

able, combined with the spectrographs of large focal ratio, which 
the light-gathering power of these telescopes makes possible, the 
loss of light at the slit becomes very serious. Thus, in the case of 
the 200-inch telsecope, if a = 1”’ (the approximate value for a star 
under good seeing conditions), F, = 30, and A = 30uyn, the ef- 
ficiency will be only 5 per cent. For poor seeing conditions or larger 
nonstellar objects the fraction of light passed by the slit becomes 
even smaller. 

From (4) it is evident that this efficiency FR is independent of the 
focal ratio of the telescope or collimator lens. Thus, reducing the 
size of the stellar image, a, by changes from the Cassegrain to New- 
tonian arrangements or by the introduction of condensing lenses in 
front of the slit produces no change in the efficiency except for 
minor differences due to reflection or absorption losses in lenses or 
mirrors. This result may also be obtained in a more fundamental 
way from thermodynamic reasoning. Thus, one may show from the 
second law that it is impossible to decrease the area on which a given 
amount of radiant energy is concentrated without correspondingly 
increasing the size of the incident cone of light; i.e., as indicated in 
(1), the focal ratio of the incident light /, must decrease proportion- 
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ately to the diameter of the image, a. Since, if the resolution is to 
be maintained, 6 must decrease proportionally to F. and therefore 
to F,, it is evident that the ratio of slit width } to image size a can- 
not be increased, and therefore the efficiency cannot be improved. 


II. DESIGN OF IMAGE-SLICER 
Since changes in area of the image are ineffective in increasing the 
efficiency, one is led to the remaining method of changing the shape 
of the image without altering its area. The optimum shape is obvi- 
ously a long, narrow band whose width is equal to that of the slit. 








If this can be accomplished by reflections at plane surfaces, the 
angle of the cone of incident light is not altered in any way, and 
therefore the equality of F, and F, is not disturbed. 

This suggests the use of a system of plane mirrors which slice up 
the image into a series of strips which are then placed end to end 
along the slit. Two views of such a system of mirrors are shown in 
Figure 1, a and b. In this diagram the beam of light is assumed to 
come from a lens of large focal ratio and, consequently, to approxi- 
mate a cylinder of uniform cross-section, equal to that of the image, 
in the immediate neighborhood of the focus. This beam, bounded 
by the rays AM and BN, is intercepted just before it reaches its 
focus by the 45° mirror MN, which is placed near one end and slight- 
ly to the side of the spectrograph slit AL (Fig. 1, b). This mirror is so 
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adjusted that the reflected beam bounded by the rays MOQSUW 
and NPRTVX is parallel to the plane of the slit jaws and just grazes 
their surface (Fig. 1, a). The mirror MN is rotated about the initial 
direction of the beam (AM or BN) into such a position that the 
reflected beam crosses over the slit KL, making a small angle @ with 
it (Fig. 1, 6). A series of mirrors, OP, OR, ST, UV, and WX, whose 
widths are equal to that of the slit, b, are placed at intervals along 
the slit where the reflected beam crosses it. They are made parallel 
to the first mirror, MN, and are therefore inclined at 45° to the 
beam. As the reflected beam passes over the slit, each of these 
mirrors intercepts a strip equal in width to the slit and reflects it 
downward through the slit in a direction parallel to the original 
beam. 

The mirrors placed over the slit have the same width as the slit 
(0.02-0.5 mm) and a length of 0.5~-2.0 mm. It is, therefore, not 
feasible to make and accurately mount mirrors of just this size. 
Fortunately, the light does not pass over the slit, and consequently 
the mirrors may be parts of surfaces that extend as far as is con- 
venient in the direction above the slit in Figure 1, 0. 

Figure 2 indicates a design for mounting these mirrors, of which 
the parts that are used (i.e., that correspond to mirrors OP, QR, 
ST, UV, and WX in Fig. 1) are indicated by the shaded areas. The 
reflecting surfaces are evaporated aluminum films applied to glass 
plates whose thicknesses are equal to the perpendicular distances be- 
tween the mirrors of Figure 1. The reflecting surfaces may therefore 
be properly spaced by stacking the glass plates in contact with each 
other. These plates must be shaped and mounted in such a way as 
to satisfy the following conditions. (1) The separate slices of the 
image, JP, PR, RT, TV, and VX, should have a length equal to the 
diameter a of the star image, and a width equal to that of the slit 0. 
(2) The successive slices should be laid down end to end on the slit 
with the sides of successive sections forming one continuous line, 
i.e., without jogs or spaces between adjacent slices. (3) The direction 
of the beams entering the slit should be parallel to the initial beam 
AM, BN. (4) No parts of the glass plates except the aluminized 
areas should interfere with either the incident or reflected beams. 
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From the foregoing conditions the following details of design can 
be deduced at once. 

1. The thickness 7 of the plates should equal a V2, in which a 
is the diameter of the star image. 

2. The number of plates should be equal to a/b in which 6 is the 
slit width. 

3. As has been already stated, the plane DPE of the mirror 
should make an angle of 45° with the plane EPK of the slit jaws. 

4. Z KPE,i.e., the angle between the slit K P and the intersection 
of the plane of the mirror with the plane of the slit jaws, should be 
1/2 — 0,in which 6 = sin 6/2 T = sin b/a. As noted above, 
6 is also the angle at which the first reflected beam crosses the slit. 


























Fic. 2 


5. The angle between the surface of the mirror DPE and the 
surface CPE should be 45° or less. J 

6. Z DPE should equal 7/2 — ¢, in which tan @ = (tan 0)/V2. 

7. The angle between the plane of the mirror DPE and the plane 
CPD should be 1/2 — ¢ or less. 

8. The plane of the mirror MN should be parallel to p!ane DPE. 

9. The focus of the telescope should be adjusted to come at the 
point in the first reflected beam MOQRST, where it strikes the 
center one of the small mirrors. 

If the cones of light reflected by the small mirrors OP, QR, etc., 
are exactly parallel, they fail to coincide at the collimator lens by an 
amount equal to the distances between these mirrors. In case the 
distance between the end mirrors becomes an appreciable fraction 
of the diameter of the collimator lens, this may result in a significant 
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loss of light. However, by making the mirror plates slightly wedge- 
shaped, with an angle a/2f. (a = image size, f. = focal length of 
collimator) between the faces EPD and GCF, the mirrors are given 
the correct tilt to bring the cones into coincidence at the collimator 
lens. In this case the center mirror, only, makes an angle of exactly 
45° with the plane of the slit jaws. 

The mirror MN is obviously not an essential part of the image- 
slicer and may be eliminated in cases where it is feasible to mount 
the spectrograph with its axis at right angles to that of the telescope. 

By following the foregoing procedure it is possible to make prac- 
tically all the light in the image pass the slit, except for the necessary 
loss at the two (or one) reflections; i.e., the efficiency of this system 
can be brought up to from 75 to go per cent, as compared with the 
5 per cent obtained in the example cited above. The foregoing 
process, however, operates by lengthening the section of the slit 
that is illuminated rather than by increasing the amount of light 
passing the unit length of the slit. The image-slicer by itself does 
not, therefore, increase the intensity of illumination on the photo- 
graphic plate, but instead gives a much wider spectrum. Thus, in 
the example of the 200-inch telescope cited above, the spectrum will 
have a width of 3 mm as normally used and of 18 mm when the 
image-slicer is introduced. 

Since the width of the spectrum obtained with the image-slicer is 
usually much greater than is necessary, it becomes advisable to 
narrow the spectrum down with a cylindrical lens and thereby ob- 
tain the corresponding increase in the intensity of illumination on 
the plate. The next two sections contain a brief discussion of the 
use of the cylindrical lenses and the limitations imposed by the 
aberrations inherent in them. 


III. USE OF CYLINDRICAL LENS TO DECREASE WIDTH OF SPECTRUM 

As normally used, the cylindrical lens is placed in front of the 
plate with the axis of its cylindrical surface parallel with the disper- 
sion of the spectrum. In this position the focal plane of the spectral 
lines is only very slightly changed by its presence. The position of 
the cylindrical lens is adjusted so that the photographic plate and 
the exit pupil of the spectrograph (either grating, prism, or camera 
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lens) are at conjugate foci of the cylindrical lens. Since the width, 
W, of the spectrum desired is almost always very small compared 
with the diameter of the exit pupil, the focal length of the cylindrical 
lens is small compared with that of the spectrograph camera, and 
the distance between plate and cylindrical lens is very nearly equal 
to the focal length of this lens. The focal length of the cylindrical 
lens is made equal to the width, W, of the spectrum desired, times 
the focal ratio FP, of the spectrograph camera lens. The-width of the 
cylindrical lens should be s'ightly greater than the width of the 
spectrum formed before its introduction. Since it is often desirable 
to reduce the width of the spectrum by a factor of from to to 20, 
cylindrical lenses of low focal ratio are required. It is, therefore, 
necessary to investigate the aberrations introduced by their presence. 


IV. EFFECT ON DEFINITION OF SPECTRAL LINES OF ABER- 
RATIONS INTRODUCED BY CYLINDRICAL LENS 

The calculation of the aberrations of a cylindrical lens depends on 
the following easily proved theorem.’ Jf any plane be drawn through 
the normal to a surface between two media, the sine of the angle between 
this plane and any ray drawn through the foot of the normal is to the 
sine of the angle between this plane and the refracted ray inversely as 
the ratio of the indices of the two media. 

Through the cylindrical lens ABCD in Figure 3 pass a plane EF 
normal to the axis of the cylinder and therefore cutting all surfaces 
of the cylindrical lens normaiiy. ‘The plane EF therefore intersects 
the photographic plate KZ in a line parallel to the spectral lines. 
Allow any ray GHTE to pass through the lens. Since the plane EF 
is normal to all surfaces, the sine of the angle between the ray and 
the plane EF must at all points be inversely proportional to 1, the 
index. Therefore, the total distance d that the ray approaches the 
plane EF is 


é l 
d=sinB D)-, (5) 


in which @ is the angle between the ray and the plane EF in air and / 
is the distance along the ray. 

' J. P. C. Southall, Principles and Methods of Geometrical Optics, p. 30, New York: 
Macmillan Co., 1910. 
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Consider two rays coming from the same point on the grating and 
making the same angle 6 with the plane EF in air. One of these rays 
passes through the center of the cylindrical lens, the other through 
the edge. The geometrical path from grating to plate, }°/, through 
the center of the lens is shorter than that through the edge of the 
lens by an amount D. 

If the surfaces of the cylindrical lenses are aplanatic surfaces, the 
optical length }°nl of the two paths must be equal. In order to satisfy 
this condition, it is necessary that the central ray pass through 
D/(n — 1) greater thickness of glass than the outer one and that 


K A B G 
Lo 
| H 














FIG. 3 


the ray through the edge pass through nD/(n — 1) greater thickness 
of air than the central one. In passing through D/(n — 1) of glass 
the central ray approaches the plane EF by the distance (sin 8) 
D/n(n — 1), while the edge ray in passing through nD/(n — 1) of 
air approaches it by the distance (sin B) nD/(m — 1). At the focus 
these rays are therefore separated by a distance sin B [nD/(n — 1) — 
D/n(n — 1)] = (sin B) D(n + 1)/n in the direction normal to EF, 
i.e., normal to the spectral line. 

If the photographic plate, and therefore the axis of the cylindrical 
surface of the lens, is normal to the axis of the cone of incident light, 
the maximum value of sin 8 is } /,, which means that the width of 
the spectral lines, 6, due to this aberration is 


_ Din+1) 
alien 2Fn * (6) 
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Case 1. Single thin cylindrical lens.—The difference in geo- 
metrical path between the central and edge rays of a single thin 
cylindrical lens of focal ratio Fx and of focal length WF, is 


_ WF, 
~ «8F 





D 


and ’ 


oo (m+ 1)We (7) 


16nk 


The minimum value of Fx that may be used, if no aberrations 
greater than 6 are permissible, is for normal incidence 


7-2 |(n+1)W 
" + \ ne (8) 





Forn = 1.5,W = 1mm,and6 = 25 uthisgives Fx = 2.04. This 
represents about the minimum value of Fx that can be obtained with 
a single cylindrical lens. If the lens surfaces are circular cylinders, 
approximately the same limit is set by the aberrations, parallel to the 
spectral lines, becoming an appreciable fraction of W at Fx = 2. 

Case 2:—Aberrations of a system of cylindrical lenses.—¥rom the 
foregoing discussion it is evident that the aberration across the 
spectral line is proportional to the difference in geometrical path 
of the rays through the center and edge of the cylindrical lens. 
Qualitatively, it is at once evident that the path difference can be 
cut down by dividing the refraction between two lenses. Thus, in 
Figure 4, if the refraction takes place at a single thin lens BB’, we 
have the path ABF for the extreme ray. On the other hand, if the 
refraction is split between the two lenses EE’ and GC’, the edge path 
becomes AEGF. This, obviously, is shorter and consequently differs 
less from the central ray CF, although the focal ratio Fx = 3 sin 
Z CFB remains the same. 

(Quantitatively, it can easily be shown that the difference in path 
D between the central or axial ray and the edge ray, which are as- 
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sumed to be originally parallel, is }°m tan y/2, where m is the dis- 
tance to which the edge ray approaches the axis and y is the angle 
between the edge ray and the central undeviated ray. Since }*m is 
the distance between the rays before the first lens is reached, it is 
evidently important to keep y as small as possible at all points. 
For a spectrum of finite width it is impossible to decrease D indefi- 
nitely by making FG’ approach zero and G’E’ become very large 
since the width of the cone of light from EE’, which increases with 
G’E’, becomes greater than the width of the lens GG’, which must 
decrease with FG’. However, by making the lens GG’ a very thick 
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plano-convex lens whose plane face is placed in contact with the 
photographic plate, is it possible to almost completely eliminate the 
path, where y has its maximum value of sin" 5 Fx, over as wide a 
spectrum as may be desired. If, following this design, the lens GG’ 
has a radius of curvature r and a thickness (1 + e)r, and the lens 
EE’ is a thin lens of focal length n»WF,/(e + 1 — ne), the minimum 
value of Fx that may be used without the spectral line width exceed- 
ing 6 is 


Fe = 2 th a wet 1-29]. 6) 


These yield the results shown in Table tr. 

Since the width of the lens GG’ must be large compared with the 
width of the spectrum formed by EF’, it is evident that r must be of 
the order of 5 mm or more if W = 1 mm. The last two rows of 
Table 1 give the numerical values of Fx for W = 1 mm,r = 5 mm, 
6 = 25 mu, and F, = 30 and for the values of m indicated. The case 
e = 1/n eliminates all aberrations along the spectral I'nes due to 




















IMAGE-SLICER 123 


are small. Case e = 1 is 


x/ 
4 


lens GG’, and those introduced by EE 
closely approximated by the use of a complete circular rod for 
lens GG’. 

If the axis of the cone of incident light is not normal to the axis 
of the cylindrical surface of the lens, but makes an angle 6 with the 
normal plane, the maximum value of sin # is approximately sin 
6 + (cos 0)/2F;. In this case the value of the aberration 6 in equa- 
tions (6) and (7) is multiplied by the factor (2F, sin 6+ cos 9). 


























TABLE 1 
e | ° 1/n I 
; I (n+1)W it n+1[ (m —1)r W I n+1f2m— r _ 
Px : _ —— |— 4 ——_] ————- + — J |— : —— +W(2-—n 
inN\ 6 1" \ 6 F ns n in\ 6 Fgn 
Focal length of WF.n 
lens EE’ WF n WF gre ~ 
2—n 
Width of spectrum 
formed by EE’ Wn 
alone ‘ Wn Wn 
2—n 
Fx(m = 3/2)....- 1.67 ¥.42 1.30 
Fx(n = 5/3) 1.55 1.26 1.06 














Likewise, the lowest permissible value of the focal ratio Fx for a given 
6 is increased by the factor of V2F, sin 6 + cos 6 in equations (8) 
and (g) and in Table r. 


V. METHODS FOR GUIDING TELESCOPE 

Since none of the light strikes the slit jaws when the image-slicer 
is in place, a change in the usual method of guiding the telescope is 
required. One procedure is to replace the first mirror MN by a right- 
angle prism to the first surface of which has been cemented, with 
Canada balsam, a thin wedge of glass with an angle of 2°—5° between 
its faces. The upper surface of the wedge is aluminized except for 
an area equal to the cross-section of the incident beam, which at that 
point is only slightly larger than the image diameter a. Owing to the 
angle between the incident beam and the normal to the upper sur- 
face of the wedge, any light reflected by this surface can be observed 
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by a microscope placed at an appropriate position beside the inci- 
dent beam. As long as the telescope is properly pointed, only 4-5 
per cent of the light is reflected at the glass air surface of the un- 
aluminized section. If the telescope drifts from its correct direc- 
tion, part of the image moves over on the aluminized area and a 
large fraction of the light is reflected. Because of the slight refraction 
of the main beam A M, BN, in passing through the wedge, a small tilt 
must be given to right-angle prism to bring the initial and final 
beams into parallelism.’ 
ASTROPHYSICAL OBSERVATORY 


CALIFORNIA INSTITUTE OF TECHNOLOGY 
June 1938 


2 A second article will give the results of performance tests, details of mounting the 
mirrors, and methods of construction of the optical parts. 




















FRAUNHOFER INTENSITIES IN THE INFRARED 
REGION dA 8800-11830 A* 


C. W. ALLEN’ 


ABSTRACT 

Equivalent breadths, measured photometrically on spectra taken with an 18-foot 
focus spectrograph at Mount Wilson Observatory, are tabulated for 188 absorption 
lines in the solar spectra between \ 8800 and \ 11830 A. 

Systematic observations of the intensities of Fraunhofer lines 
have already been published by the author for the visual? and the 
near infrared spectrum. Since for many solar problems it is desir- 
able that such information should cover as wide a wave-length range 
as possible, the observations have been extended as far into the 
infrared as the available plates would allow. 

Photometry.—Spectra of the center of the sun’s disk were obtained 
with an 8-inch Michelson plane grating mounted at the 18-foot focus 
of the spectrograph pit of the Snow telescope of the Mount Wilson 
Observatory. The plates were Eastman 144P, 144M, 144Q, 1Z, 
backed and ammoniated, and developed for 5 minutes at 65° F with 
a hydroquinone developer of strong contrast (X-ray). Usually three 
spectra of varying exposure times were recorded on each plate, the 
best two being selected for photometry. The plates were standard- 
ized with solar spectra taken through a step-slit and with similar 
exposure times. For this purpose the intensity of the sun’s image 
was decreased by inserting a number of black cheesecloth screens 
into the telescope beam. For each measurable line the uncorrected 
central intensity, 7», and the measured breadth of the line at half 
its maximum depth, A,,, were obtained. For the majority of lines 
the equivalent breadth is then given by W = 1.22 A,(1 — fm), where 
the numerical factor is obtained from considerations of scattered 


* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 594 

' Hackett Research Student. 

2 Mem. Commonwealth Solar Obs. Canberra, 1, No. 5, 1934. 


3 M.N., 96, 843, 1930. 
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light and line contour. For the strongest lines the factor is some- 
what greater. This procedure enables allowance to be made for the 
appreciable difference in the breadth of Fraunhofer lines in the infra- 
red. The effect of this difference was not realized when my last in- 
tensity list? was published, and hence the equivalent breadths of a 
few of the broadest lines given in that list may be too small. 

An attempt has been made to correct 7, for the blending of 
atmospheric or solar lines. When it was not possible to measure A,,, 
a mean value for the region was adopted. 

With the exception of a few gaps, four plates covered each region 
from \ 8800 to A 11000 A. Beyond 11000 A, however, only one bad- 
ly fogged plate was available, and the results must be regarded as 
tentative. Where different type of plates overlapped, there was some 
evidence that plates of higher contrast gave systematically smaller 
equivalent breadths. It is hoped to investigate this question later, 
but for the present the mean is taken. The present observations are 
in reasonable agreement with those of Minnaert and Bannier.4 The 
comparison was made, however, in a region where there is no blend- 
ing with atmospheric lines, and hence the agreement may give an 
exaggerated impression of the general accuracy. 

Results.—The results are given in Table 1. The first four columns 
contain material kindly supplied to me in advance of publication by 
Mr. H. D. Babcock and Dr. Charlotte E. Moore. Some of the wave 
lengths, especially beyond 11000 A, may be further improved by 
small changes. The intensities, by Miss Moore, are not in every case 
final. They are designed to be concordant with Rowland’s intensities 
in the region near \ 6200. The fifth column contains the equivalent 
breadths, W, resulting from the present investigation; the sixth, 
the corresponding results of Minnaert and Bannier4 (MB), Dahme* 
(Da), and Minnaert and Genard® (MG). All values are in mil- 
liangstroms. The last column gives the quantity W/X, which is often 
more convenient in application than W itself.7 These values of W 
are averages from the two preceding columns. More doubtful values 
of W are in parentheses. 

4Zs.f. Ap., 11, 392, 1936. 


5 [bid., 11, 93, 1935. 
® Ibid., 10, 377, 1935. 7D. H. Menzel, Ap. J., 84, 462, 1936. 
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TABLE 1 

















W IN 

WavE MILLIANGSTROMS 
Seaien ELEMENT E.P. INTENSITY 10°W /d 

Author Others 

MG 
8804.637.... Fe 2.269 3 CO Caer enn Pe 6.4 
8806.775.... Mg 4.327 12 495 512 iy ae 
S808 1735 .«:. } | i eres c I OU eee 4.0 
8809 .406.... Ni 3.881 —2 COR PR er es) 3.6 
8824.236.... Fe 2.188 9 Be Cae. Soe 20.5 
8838 .441.... Fe 2.846 6 FOP eos shea 11.6 
8846.750.... Be  ‘Pentece. 3 CoO Eerie ie 5.9 
8862.054.... CS ibs :ecree — 2 BO Bracket erat 2.3 
8862.563.... Ni 4.072 3 Ce) PCr 6.9 
8866 .943.... Fe 4.529 8 Beilin, Poze sah. carer 17.4 
8868 .443.... | ne ner 3 CO an] Beeson 5.7 
8876.031.... Re ‘Boxcar. cet I TE Ne diewvtosee 3.8 
8883.68..... Si 5.925 | = 2 Ue ag Sewers ek: 1.5 
8884.24..... eC t«é«C | — 3 JOO Eee es ¥.3 
8892.738.... Si 5-958 4 We Nive aaieo wre 8.7 
$800: 222... Atme 8 — Jaisescse 2 AGS ViscxeKnees 5.2 
8905.989.... S ne Ocererne ° es Werte chaets 2.6 
Sor? 102... .. 2 ee, Pers 7 ) <a Spee rgrer 4227 
8920.036.... Fe, Atm? (V)*| 5.042 3 Gi be cee ted 6.8 
$027.-'570: . «. ae Ree 3 11 An) ES eee 6.8 
$024. 0052.<: Ge § + Wests — 2 Le ed Pree. a4 
8027.302. ... Gos, .... Henaeiecs 4 oe aM Ree eee 434-7 
8943.058.... Fe (Na??)* asians aes I Ath les cca ae 
8945.198.... Fe 5.011 3 oa eerie 8.9 
8947.197.... Ci,  “bovtenaee — 2 0G) Sel ear are 1.9 
$650.297:... «.. | Se Soper — 3 Oe Seer ey 
8067). 72. .... OMe — CPissimeca — 2 CIS) Vena saree 2.0 
8968.20..... | ne cree — 2 P 7 gh Pere 2.7 
8075.43... Fe 2.977| — 1 SR ira stncaes 9.0 
8976.88..... 5; Sn ae re ee — 2 SOP Pesncceeweon 2.0 
8070. 23. ..-- ee ee Loreen — 3 Se Pres yp nectins 2.0 
8984.898.... Fe 5.078 | — 1 Bae Ds tigre se 3.6 
8999.580.... Fe 2.819 2 FAG Wee ce sccinds 14.0 
Q00S:52....... |: | nn Serre ° GE) Pic aioarecces 7.6 
goog .835.... Cr 3.307 ° GR sn See os 6.9 
GOION673. «3. Be cra rats ° kel eee 3.9 
QOFA.. 60. 5... BRE 828 =6 PE ee — I tall) Cee: Wey 2.4 
Q01G: 77...6<. Fe 5.078 | — 2 Abs — Pitsid shoal seers 2.3 
OOIE. Sos ces Cr 3.309 I Ger fev. vasa 6.4 
g024.38..... me wy | OMe ° CGB loss. xa letra 6.9 


























* Parentheses indicate masked. 
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TABLE 1—Continued 

















W IN 
ae MILLIANGSTROMS 
testi ELEMENT E.P. INTENSITY 10°W 
LENGTH 
Author Others 

O020.97S.... .<. Be ‘fe chaseds — 2 21 2 
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TABLE 1—Continued 
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TABLE 1—Continued 
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TABLE 1—Continued 
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A study of intensities within multiplets and groups of multiplets 
is being made from the infrared equivalent breadths now available. 
For a complete understanding of such problems, however, it is very 
desirable that intensity observations in the AA 3000-4000 A region 
of the solar spectrum should be made, and it is to be hoped that an 
observatory with suitable equipment will undertake the study of 
this region. 
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THE THEORY OF CYCLICAL TRANSITIONS 
L. G. HENYEY 


ABSTRACT 

The general theory of cyclical transitions in the absence of collisions is discussed, 
first for a finite number of discrete states. It is found that the rigorous treatment of the 
radiation field can be carried conveniently to a fairly advanced state, short of the com- 
plete solution of the problem. It is found that there exist a number of first integrals of 
the differential equations, some of which are generalizations of the constancy of the 
net flux and of the linear dependence of the radiation pressure on the optical depth in 
the simple scattering problem. The remaining integrals are simple relationships among 
the optical depths in the various frequencies. 

The second part of the discussion is devoted to a study of the properties of the radia- 
tion field when the exciting radiation is highly diluted. The approximation appropriate 
to a high dilution is introduced, and it is found that the problem can be completely 
solved to within the Milne-Eddington approximation to the directional dependence of 
the specific intensity. The boundary conditions for a planetary nebula are introduced, 
and the values of the emergent fluxes are obtained. Physically the approximation for a 
high dilution is equivalent to neglecting the transitions corresponding to stimulated 
emissions and those corresponding to absorption when the lower state is not the ground 
state. However, these facts do not constitute the basis for the approximation; it is in- 
troduced from more fundamental considerations. 

The last part of the theoretical discussion pertains to the generalization of the theory 
to an infinite number of discrete states together with a continuum. This formal exten- 
sion follows in a natural way and depends on no new physical concepts. 

Finally there is a section devoted to some computations on the four-state problem 
for a planetary nebula. One interesting result is the weakness of the (1, 3) radiation 
(Ly 8 for hydrogen) when the fourth state is identified with the continuum. 


INTRODUCTION 


The main features of the problem of the excitation of emission 
nebulae have been ascertained through the efforts of various in- 
vestigators. Zanstra’s ideas have been widely accepted and are gen- 
erally regarded as the foundations upon which theoretical investi- 
gations must rest. The analytical application of his ideas to a study 
of the radiative equilibrium of a planetary nebula was first made by 
Ambarzumian,’ who considered what can properly be called the 
three-state problem. The problem was later considered by Chan- 
drasekhar,? who traced the conversion of the incident ultraviolet 
radiation into radiation in Lyman a through the nebular layers. 

The problem of the Balmer decrement has been studied theoreti- 
cally by Cillié’ and others. The result of these studies has been 


* M.N., 93, 50, 1931. Also Pulk. Obs. Bull., No. 13. 
2Zs.f. Ap., 9, 266, 1935. 3M.N., 96, 771, 1936. 
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rather disappointing, in that the agreement with observations was 
found to be rather poor. It had been found that the lower members 
of the Balmer series are relatively much more intense than can be 
expected on the basis of the electron-capture mechanism. An at- 
tempt by Berman to remove the discrepancy by correcting for inter- 
stellar reddening has shown that the larger part of the differences 
between observational results and theoretical predictions can be re- 
moved but that there still remains some discrepancy. 

It is the purpose of the present investigation to provide a more 
rigorous and systematic treatment of the problem of the radiative 
equilibrium involving cyclical transitions than has yet been attempt- 
ed. The discussion has been limited to a nebula in static equilibrium. 
Although Zanstra‘ has shown that such an object cannot remain in a 
static state for any length of time, it is nevertheless desirable to ob- 
tain as complete a picture as is possible regarding the state of affairs 
in such a nebula before introducing refinements. This treatment pro- 
vides the basis for future extensions and gives the correct method of 
procedure in formulating the problem. It shows, for instance, that 
the complete treatment of the continuum, or the equilibrium of the 
electrons, provides no difficulties which differ from what one meets 
in the treatment of the discrete states. 

The investigation is roughly divided into two main sections. In 
the first section we will study the general relationships concerning 
the condition of radiative equilibrium and the derivation of several 
integrals describing the field of radiation. In the second section we 
will make use of the high dilution of nebular radiations to solve for- 
mally the complete representation of the radiation field. It will be 
convenient at certain stages to regard the continuum as a discrete 
state, but eventually we will consider the correct method of handling 
this part of the problem. 

This investigation originated from a suggestion made to me by 
Dr. Chandrasekhar regarding the desirability of a complete and 
rigorous treatment of the three-state problem in the absence of col- 
lision effects. In working on this suggestion, I found that all the ex- 
plicit relationships which I could derive could be generalized to any 
number of states, showing that a start can already be made on the 


4 [bid., 97, 37, 1936. 
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more general n-state problem; and, although the work has passed 
beyond the original suggestion, I am still greatly indebted to Dr. 
Chandrasekhar for it. 
THE EQUATION OF TRANSFER 

Let A;;, B;;, and B;; be the Einstein probability co-efficients of 
spontaneous emission, stimulated emission, and absorption, re- 
spectively, defined with respect to the intensity of the radiation. 
If 1; is the number of atoms in the state (j) per unit volume, the 
number of downward transitions from state (7) to state (7) is 


I ~ 
n(Ajit pa B;i{ I, dw) (1) 


per unit time and unit volume. The number of upward transitions, 
under similar circumstances, is 


I 
an n; Bijf Iv dw . (2) 


We suppose that the emission and absorption take place within a 
frequency interval Av;;, and we consider the average emission and 
absorption coefficients. The average atomic emission coefficient is 


hv 
4m Apv;; me (3) 


the average atomic stimulated emission coefficient is 


hv 
4m Apv;; Bit (4) 


and, finally, the average atomic absorption coefficient is 


hv 
4m Ap;; me (5) 


Hence, the equation of transfer becomes 


11, h 
7 be poe (njAji + njByl, — n Bil) . (6) 
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We can now make use of the familiar relationships between the vari- 
ous Einstein coefficients, namely, 


Bi? gj 
Aji = 2hv3 gi : (7) 
B;; os Gc? 

A ji a 2hv3 : (8) 


The quantities g; are the statistical weights. Eliminating the B’s 
from (6), we find that 


dl, hy nN; nm; n;| ¢? 
ds ~ 814i 4m Av ;; & i fe gi 2hvs 1) 


In the discussion which follows, it is more convenient to describe 
the radiation field in terms of the flow of photons rather than in 
terms of the energy flux. Furthermore, we shall use the total inten- 
sity over the entire line rather than the intensity per unit frequency 
interval. We therefore introduce the symbol J;; defined by the 
equation 


(9) 


where the subscripts refer to the lower and upper states, respectively. 
As will become apparent, the ratios ;/g; play a more fundamental 
role than the n’s themselves. We denote this ratio by Nj, and we will 
refer to it as the number of atoms in state (7) per unit volume per 
unit weight. Finally, we make the abbreviations 


I 
ra giAji = Bij, (10) 
Cc? gi Aji = i 

Srv? Avi; sl ) 


After all these changes the equation of transfer finally becomes 


d dij e i : 
oe BiiN; + (Nj — Ni)picdii (12) 


for transitions between discrete states. 
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In the case of the bound-free or free-free transitions it is possible 
to write the equations of transfer in an analogous form. In order to 
clarify the analysis, we now disregard the continuum and assume 
temporarily that an atom has only discrete states, and for further 
simplicity a finite number—say r. Later we shall see that the gener- 
alization of the theory offers no difficulties; and, when we do con- 
sider these refinements, we shall see that the results follow directly. 


























Fis Je 
1.0 
. aA 
0.0 
° I 2 3T 
l'ic. 1.—From top to bottom the curves correspond to py: = 1, 3, 3, }. These 


curves also represent (1 — py) GaP, 2 — Pa) (Py2) Fag, and (1 — py) (P43) 2° 


THE CONDITION OF RADIATIVE EQUILIBRIUM 


In order that we may solve (12), it is necessary that we evaluate 
the N’s. Now, the population of the various atomic states is deter- 
mined in a steady state by the condition of radiative equilibrium. 
In the case when the gas is mechanically in static equilibrium, this 
condition is simply the mathematical expression of the fact that at 
any point the number of transitions into a state during an interval 
of time is equal to the number of transitions out of this state in the 
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same interval of time, all transitions taking place spontaneously 


or under the influence of radiation. If we write 


-, dd 
Jii = fdii = ’ (13) 


where the integration is performed over all directions, we have for 
the condition of radiative equilibrium 


SS Nimiigir a im N (Bis + wisi) 
j= j=it+i 
= SN Bii + Hii fii) + ba N mij! ij ¢ 
jz i=i+1 


Rearranging the terms, we have 


2. (N.Bii + {Ni - Nj} mii fii) | 
j=t ; | (14) 
- a. (NB:; + {N; eo Ni} miifis) =O 
j=itt } 


The equations in (14) are not all linearly independent, for, if we 
sum over all possible values of (7), we get identically zero. The equa- 
tions serve in principle to calculate the ratios of the N’s in terms of 
the intensities in the radiation field. Obviously, an explicit solution 
of (14) will be exceedingly complicated. However, various general 
integrals of the differential equations of the radiation field can be 
obtained without solving the equations of radiative equilibrium. 

The operator d/ds, occurring in the equation of transfer, can be 
expanded into the well-known form 


” 


0 0 0 
l - m. n—, 
Ox * oy sl Oz 


where (1, m, 2) are the direction cosines of the ray. The vectorial net 
photon flux 7Y;; has components given by 


1 1 cs a 1 
~ Teas Jil = ) Fy, jf Si m ’ Fe, Go= f Dijn = . (15) 
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An integration over all solid angles of both sides of (12) leads to the 


result 
r div Fij —_ ( ) 
4\ Ox dy Oz (16) 
= BiiN; a (N; ion N i)miifii . 
Consequently, by (14), 
div (Sis: -»> 5, =o. (17) 
\jJ=1 jJ=t+I / 


If the material is stratified in parallel planes, (17) becomes, for the 
fluxes normal to the plane of stratification, 


dz\ <_< — 
j=l j=iti 


: (Ss: “* Ss.) =o, 


where z is the co-ordinate whose variation is normal to the planes of 
stratification. Hence, 


i—! r 

~ ~ 
> Fi - > Fi: = Ets (18) 
j=l j=iti1 


If the system has spherical symmetry, we can deduce from (17) that 
for the radial fluxes, 


~ ™~ a; 
DF 2 Fa = 5s (19) 
j=tI j=iti 


where 7 is the distance from the center of symmetry. In general, 


iI r 


~ ~ 
DF - DY Fi = curl As, (20) 
j=! j=itiI 


where A; is an arbitrary vector field. These relationships, which we 
will call the flux integrals, constitute the generalization of the 
familiar constancy of the net flux in the pure scattering problem. 
Not all of these integrals are independent, for, if we sum them over 
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all the states, we have identically zero. Hence, there are r — 1 in- 
dependent relationships. 

While the equations (18), (19), and (20) are concerned with the 
photon flow, they guarantee the conservation of energy in the radia- 


Fi; Bf’, X 105 
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FIG. 2.—py, = }. From top to bottom we have 4. and p;: equal to (0, 3), (4, 3), 


(0, 3), G, 4). 


tion field as in the simple scattering problem. The total energy flux 
is proportional to 


aim. r 


a | 
\ > Vij ij . 
i=1 j=i+1 


From the Bohr frequency relation this equals 


F=3 =f F=f 
™~ "Sie | ~ 
>. > (viy — Vir) Fij + Do vieFir . 


i=1 j=iti s=1 
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Rearranging the terms in the double summation, this reduces to 
~ / i-! 
‘ 
> n(- > Fiji + Ss 5, . 
i=I ry 1 j= =. / 


From (17) it can be-seen that the divergence of this expression is 
zero. This proves the proposition. 
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FIG. 3.—py = 3. From top to bottom we have p,2. and p;: equal to (0, 3), (4, 3), 


(0, 3), (i, ; 


The optical depth 7;; along a ray in the (7j) radiation is defined 
by 
dri; = (N; = Ni) pijds . (21) 


For any three states (7, 7, and k) such that i <7 < , it follows i im- 
mediately from (21) that 


a(7 st Tik ra) — 
bij Mik Mik 
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v4 


along the ray. Or, 


Tij + T ik _ Tik — ( 
8 ’ 
Kij Mik Mik 


No 
to 
wa 


provided they are all measured from the same point. We shall refer 
to these as the r-inlegrals. Now, if we write (22) for all possible 
sets of three states, not all the equations will be independent. In 
fact, there are only [(r — 1) (ry — 2)|/2 independent integrals. To 
see this, let us consider the set of equations 


T1j T jk Tih 


——=0. ics te Se). 
Mij Mik Mik 

This set contains only such equations as are independent, since only 

one of them has the term 7;/u;, and it contains exactly [(r — 1) 

(yr — 2)|/2 equations. Now, any combination of three states other 

than those represented in this set leads to an integral which can be 

derived from those in the set, for 


Taj . Tjk _ Tik a (2 re Tig =i) fe 6 -- Tik =) 
Kij Mik Mik Mii Ki} Mi Mij Mik Mik 
_ (7 Tik — Tik 


2 QQE<ci<j<k <4), 
\Mii Mik Mik 


proving the proposition. 

Now, suppose that J;; represents the intensity of the incident 
illumination after passage through a portion of the nebula. Then, 
along a ray, 


dd; ; = 
ds = (N a N i) Midi . (23) 
Hence, 
By = HM, (24) 


where the zero superscripts refer to the surface at which the incident 
radiation enters the nebula. The 7-integrals can be written in the 
alternate form: 

Ji AF k IT Mii: 


I j I 
ma - + — ee ( 
nis 8 Je, + win 8 Jp, woe 8 Thy 


No 
un 
~~ 
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In a preceding paragraph we considered the generalization of the 
flux integral in the simple scattering problem. By virtue of these 
relationships we are able to derive one other integral which is a gener- 
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alization of another relationship in the theory of simple scattering— 

the well-known’ linear dependence of the pressure on the optical 

depth. The possibility of generalizing this integral was suggested to 
5 Cf. Milne, Handbuch d. Ap., 3, Eq. 125. 
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the author by Dr. Chandrasekhar, who found a special form of 
the generalization in his study of the three-state problem. In deriv- 
ing the integral we limit ourselves to the case of plane-parallel strati- 


fication. Then we have 





add; ; 
cos @ ie = BN; + (N; = Na) wij dij ° (26) 
Since 
I si : 
Jii — 2 { Ji; sin dé P 
$3 = 2 Ji; sin 6 cos 6dé , (27) 
Ki; = f Qi; sin 6 cos? 6dé , 
we have 
dF i; 
1 asi , - ; 
— = BiN; + (Ni — Niwifr, | 
‘= | (28) 
dK i; Aare 7 - 
dz a 4(N; on IN i) mij F ij . 


From the second of these, 


r r r r 
qdxyuvyii : = =. ‘ 
N N\ Ki; = \ \ (N; — Ni)Fi;. 
dz od fod hij 4 = Reurared 
i=1 J=it+I1 i=1 J=1+1 


After some reduction the right-hand side becomes 


a a 
: oN > Fi - >) Fi 


- j= j=tti 


) 
| 
\ 
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From the flux integrals this equals 


r ( i-t r 


I ae = ae ~ 
SN SH, -— SHS, 
i=I j=!1 j=iti 


where the superscripts refer to the values of the quantities at some 
arbitrary but fixed value of z. If we reverse the reduction, this be- 
comes 
roof 
; S Swi - v)H,;- 


t=1 j=i+1 
From the definition of the optical depth the equation becomes 


r r r if 
d a ne | 1d SQ Ti 

* Fx, - > + -H.. 
dz met Uj; 4d mt med Uj 

#=1 j=i+1 i=1t j=i+r 


Integrating, we finally have the K-zn/egral, 


~ & 
‘\ (Ki; — $92;7:;) = constant . (29) 
ad Mij 
1,J 
In concluding this section, we remark that we have found r — 1 
independent flux integrals, [(r — 1) (r — 2)]/2 independent 7-in- 
tegrals, and one K-integral, giving a total of (7? — 3r + 4)/2 inde- 
pendent integrals. 
THE THREE-STATE PROBLEM 
As an example of the preceding analysis and for use in later sec- 
tions of this paper, in this section we formulate the three-state prob- 
lem for a plane-parallel model with perpendicularly incident radia- 


tion. 
The condition of radiative equilibrium gives rise to the two equa- 


tions: 
(iia Pro + Hrs us) 7 (B:.2 + His fr2) cs (B:; + Hiss) en Oks 
es MifoN, + (B.. ++ Mra /pi2 + 23 a3)N2 Eom (B25 - M2323) N 3 =o. 
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Solving for the ratios of the V’s, we have 


N,« B..(B:; + B.;) + (B:; + B23) Mr2/pr2 
+ Br. bis fis + (B:. + B13) Mas fas + Mrz brs Pra fas 

+ Hizb; Pra fos + Hr3 Mas P13 $23 ’ | 

\ 


N, a (B:; + B23) Mi2/pr2 + B23 Mrs rs + Mrz brs Profs | 
+ M12 M23 Jr2 fps + M13 Mas Pus P25 ’ 


N, oa Br2 i313 + Mrz Mrs Pras + Miz Mas $r2 fas + M13 Mas P1323 . 


(30) 


Next, we need the quantities 


N,—N,« B..(B:; + B.;) + (B:. = B23) x3 $13 | 
+ (B12 + Bi 3) Mas (p23 ) | 


N, ack N; a B,.(B:; + B.;) a (B:; + B:3)bi2ffr2 (31) 
+ (Bi + B53) Mas fpos ’ | 


N,- N; a (B:; + B23) Mr2fr2 = (B,. = B23) ors rs , 


| 


From these equations and (28) 


I dF 2 I dF; I dF, 
4 7” =— ‘ rT = 3 a oc BB 25 bis /p13 + B25 bi2bis Profs 


7 Br 3 Mrz Mos Pra fpos + Br2 brs Mas $rs 2s | 


(32) 


and 


dK: 


dz 


i = 1! B.2(B:; + B.;) + (B,.2 a B.;) ois fas | 
+ (Bi. + B13) Mas Jos} Biot i2 > | 


dK; 


dz 


ae , {B:2(B,; + B.;) + (B:; + B23) Mr2/fr2 (33) 
+ (B:: i Br 3) bos 2s} Mise 13) 
dK:, 


~ 


Sy 41 (B:; + B23) Mia /pr2 — (@: — 1B 25) Mrs (Jr3 } Mas Fas . 
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The normal optical depths are defined by 


1 dt ao — B:.(B:; + B.;) — (@. — Bos) M33 | 


Mi2 dz | 
ome (Bx. + Br;) Mas fos ’ 

I rs aco — B..(B:; + B.;) ~~ (B:, + B.;)Mi2fr2 (34) 

Mi3 GS | 
= (B:. + Bi 3) Mas fas ; | 

dr», | | 
I = & <- (B:; + B.;)bi2$r2 ~ (B,. a B23) Mrs $13 ‘ | 
M23 @3 


The factor of proportionality is the same in all of the preceding 
equations. Equations (30) show that this factor is always positive 
and different from zero, so that, if we denote it by A, 


du 


dz 


determines a function #(z) which is monotone increasing. Hence, we 
can use it as the independent variable; and by substituting it for 
z, we can replace all of the proportionality signs in equations (32) 
to (34) by equality signs. 

In order to complete the formulation of the problem we must in- 
troduce the Milne-Eddington approximation. If there is incident 
radiation which is concentrated in a small solid angle, we must sepa- 
rate it from the diffuse radiation before we apply the approximation. 
If we denote the incident photon flux by 7d;;, and if the rays are 
along the positive z direction, we have that the J;; and the K;; for 
the incident light is equal to 3;;. Therefore, the correct approxima- 
tion is 


3(Kij = 423i) = Jii a iDii » 


or 


Ki; cae api + §Dij ° (35) 


In general, when the third state is identified with the continuum, 
the diffuse radiation will be determined largely by 3,,. It is of inter- 
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est to consider the case of both 3,, and 3,,; equal to zero in the plane- 
tary problem, since we can establish another integral under these 
circumstances. From the flux integrals or from (42) we see that &,, 
and &,, differ only by a constant. Since the diffuse fluxes must be 
zero on the inner boundary of a planetary nebula, it follows that this 
constant difference must be zero. Therefore, by straightforward 
substitution we have from (32) and (33) 


1 (B B F dF. 3 ( IK;: B ( 
16 (; “ees as) 12 du + a 23 fe da 7 A Jo + 


i ec 
Mo;/ du 
Using the approximation (35), when 3,, = 3.,; = 0, we have 


\ B.;\ 2 
eS . p.,)F * #. fir — BAP a Mas ) 


= constant . (36) 


A special form of this integral has been derived by Chandrasekhar. 


THE NEBULAR PROBLEM 
In the preceding discussion we have not concerned ourselves with 
the essential property of nebular radiation, namely, its high dilution. 
If the incident light is roughly distributed according to dilute black- 
body radiation, this high dilution implies that 


2 


hy | ae a 


except in the extremely long wave-length region of the spectrum. In 
our notation this inequality becomes 


w 
~I 
—S 


Bo <1. (. 


Now, in the nebula the incident radiation is redistributed over differ- 
ent parts of the spectrum, the most striking being the enormous gain 
of Ly a in the diffuse radiation over the direct. This gain probably 
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does not exceed to? or, if the Doppler broadening of the lines is 
lower than is usually accepted, 10°. However, this gain is hardly 
sufficient in the great majority of planetary or diffuse galactic nebu- 
lae to neutralize the smallness of the dilution factor, which lies in 
the range <10°.° Even in Be or Wolf-Rayet stars, where the dilu- 
tion factor is not as extreme asin nebulae, an important first approxi- 
mation to the actual situation can be worked out on the basis of 
(37). The importance of the inequality lies in the fact that it enables 
us to solve the conditions of radiative equilibrium explicitly in a 
very satisfactory form. 

An examination of the complete solution (30) for the three-state 
problem gives the clue to the form of the general solution. Applying 
(37) and remembering that B,, ~ B,., as far as order of magnitude 
is concerned, we have 


B. 23 M13 
B., +B) 2? 





N, = 3 
N,= oe +B 





a) Fon Br y) 
N, By + B37" " 


These values suggest the following generalization (j > 1) for r 
states: 
] I _ Mr; ty . Mil 8 
N, = 5 B,,o" + 2) Vit g du: (3 ) 
l=j+1 
To determine the coefficients in (38) we substitute these equations 
into (14). Retaining only first-order terms, (14) (for 7 > 1) be- 
comes, after some rearranging, 


( i~! ) 
| ie ek 
Benth 
j=t ) 
i—1! l—1 ) 
~ | Vil ~ Biyi ~ Biivin 
mS Miia 4 ) \ Bi . i SN } =O 
B:. —— B: —_ B:. | 
j= ea j=! j=it+i ) 


6 This argument is essentially due to Ambarzumian, M.N., 95, 469, 1935. 
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Each of these equations is satisfied by equating the coefficients of 
each of the ¥’s to zero. Doing this, we obtain the following recur- 
rence formulae for the coefficients in (38): 


B:; 
i = ) (39) 
S*B;: 
— 
jut 
( l—t ) 
I ' 
WAL = ae Biv + > Bisvin f . (40) 
eS { j=i 
Bi: ’ Jj +e ) 
j=l 
Now let 
Bi; B:; 
Pii — soe ’ Pi: ae aes = Vis (41) 
j j 
' — 
> hy Bi; 
k=1 k=1 
and let 
B:1 
va. * qui - (42) 
‘2. 
a 
k=! 
From (41), 
j-1 
Se 
aes 


i=I 


We can readily see that /;; is the probability that an atom, under- 
going a spontaneous transition from a state (/), will pass directly into 
the state (7). From (40), (41), and (42) 

j-I 


qui = put > ij Pii - (43) 


j=iti 


Now consider the probability, analogous to p;;, that the atom, after 
leaving (7), will arrive eventually in (7) after an arbitrary number of 
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downward transitions. But these probabilities must satisfy recur- 
rence formulae identical with (43) if we identify them with the q’s. 
Therefore, g;; is, in fact, this quantity. 

Equation (38) can be written in the form 


r 


Nie . 
SB: = Mai fps + > qijMri fpr . (44) 


N 
+ p= l=j+1 


Equation (44) is precisely what we should expect if the effect of 
transitions, corresponding to absorption and stimulated emissions 
and not involving the ground state, is negligible. For the terms on 
the right-hand side give the frequency with which atoms arrive into 
state (7), owing to absorptions in the principal series, whereas the 
term on the left represents the frequency of departure. This ex- 
planation of (44) is the underlying principle in recent calculations of 
the theoretical Balmer decrement. 

Having determined the population of each state, we can now cal- 
culate the intensities in the radiation field. The first step is the in- 
vestigation of the transitions involving the ground state. The emis- 
sion coefficient for the case of a discrete upper state is 


{ 


{ “e ] 
B,;N; = N; pin § Mairi + > qiibri fui . (45) 


l=j+1 
The absorption coefficients are 
(Ni — Nj)ini ~ Nui; ; 


by virtue of the fact that V;/N, is small. Hence, for a medium 
stratified in parallel planes, 


I dF; 


= p 
4 dl Le SS Bim ge 7. Miffsi » 


k=j 


where 
t=. JNidz 
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By; = QkiPir 5 Bij = Pir - (46) 


It is convenient to separate the incident and the diffuse intensities. 
Let us denote the quantities that refer to the total intensity by 
means of an asterisk. Hence, 


\ 


Pi = Ju +3du, | 
Fi, = Fij + Bij ’ } (47) 
Ki, — Ki; + 1 Biss 

where 7d;; is the net flux in the incident radiation. We have for 


our complete system, using the Milne-Eddington approximation 
and dropping the subscript 1, 


I dF; - ~ 1g ( 
4 dt = ae Poste tT k) = wif ) 

d ; | > 
d; ae tMiF ) (48) 
at 
d3; j 
iC 


Since the equations are linear, we can make another simplification. 
We suppose that all the incident intensities are zero except 3,. 
The complete solution is simply the sum of all these special solutions 
with 7 taking all possible values. 

Eliminating ¥; from (48), and writing 


N} = 3K; (1 9 B;i) ) 


we have 


d?9; : 
d: Nii wa > Bij bic rE 1 Brj Mj Mrdr : (49) 


dt? ' 
k=j+1 
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We try a solution of the 


form 


3H . S Qj 3H 
gi = A008 + Bio + ws at Harb (50) 


d; 


’ 


k=j+1 


Now 


k “~ 
— Ne = (M — A) Je — 3me ~ Biro fn — fukBridy « 


l=k+1 


wy 
t 


dt 


Using this equation and introducing (50) into (49), we find that 


k—t1 


d’y ue" | ~~ 

Oo! 2( a N a. Qj} 2 No XN \ 

dt? dj Ji ee é Mi uz (Xj d; ) ¥ . aj Bx 
P= ji =j+1 


/ 
; ad; 


2 
. > > _ 
+ TMiMrdd: ee (us; — d;) =o \ Ayr Brk t . 


‘ k=j+1 
To satisfy (49) we equate the coefficients of J, and 3, in (49) to those 
in the preceding equation. From this we derive the recurrence for- 


mula 
F k—-1 
3uc OC i aa 
Li. = — 1 b. 5 
Qik ue he By; 3 Qj Bur | (51) 
: l=j+1 
and the formula 
, r 
br a 
= , — 1 r x 52 
a; if — pe) - aj. Bri (5 ) 
br i 


It follows that the solution of the differential equations (48) can 
be expanded into (50) provided that no two of the )’s are equal 
and provided that no \ is equal to u,. When these conditions are 
not satisfied, it is necessary to introduce polynomial terms in ¢ into 
the solution. We will not consider these degenerate cases separately, 
since, as is well known, their solutions can be derived from the non- 
degenerate solution after the boundary conditions have been as- 
signed. 
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By means of the second equation in (48) we can derive the fluxes. 


We have 


r 
ae 
Fj = YX anFe + a8, — gf Ape — Be}. (53) 
k=j+1 


We have now to introduce the boundary conditions. Let ¢, be the 
value of ¢ at the outer boundary of the nebula. Then we have the 
well-known conditions for a planetary nebula: 


oO, 


ty. 


J; =0 when t 


Fj = 29; when t 


F] 


The first condition gives 


a; “ 
Aj — By=" Be, 


while the second gives 


I ; } j 
Ss ai( ¥ ~ 1) Fa(t) + ta,{ 3 - 2) Br(ts) 
2 a Mk r 
k=j+i 


. Shi 3H) 
+ A,( r, a 2) 4+ B,( d, al 2)e Ajé =oOoO. 


Solving for the constants, we find that 








aye (9 - 2 emt —2 N ae? — t ) F(t) _ a,( 3H _ 2)3,(t) 
I dj Ps od Kk Mr 
eas JtI - 
a 7 — H \ 
(Bi sone (= adem 
a5 ( 5H + 2)ent + 2 v ae ( M! = is )Fu(t) + a,( 3 = 2) S,(t) 
I i Phy v's Bk \ Mr 
B; = — AS ee icant 
: (B+ sone (Bi a)em 

















ee 
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The emergent flux becomes, after substituting the values of the con- 





stants, 
3 2 ee ) 
r . cosh Ajt; + ~— sinh Ajt, 
~ 4 > 
F (ts) = YS ai Fi(t) i = 
k=jt+t 3 cosh Ajt, + sinh A jt, 
\j Mj ( ) 
. , ’ ( (54 
cosh Ajt, + — sinh Ajt, — “= ert: | 
+ ee 2 ace : 
> cosh Ajt; + ~ sinh Ajt, 
Mj 





When j = 2, B.. = px = 1 and X, = o. From (43), (51), and (52) 
ada = I, GS —~ 1s 


Equation (54) becomes 


~ 


F(t) =. = Fie(ts) T SP (1 cs erty) ’ 


M@- 


= 


=2 


] 


which also follows from the flux integral. 

We are now in a position to determine the intensities for transi- 
tions not involving the ground state. The emission coefficient B;;N; 
is, by (46), 

r 
N, i BijiMs (Yr + idk), 


k=) 


where we have written 


Brit = QeiPji = Bei * (55) 


The absorption coefficient 
ui(N; — Nj) 


is a linear combination of intensities and is consequently a small 
quantity. When it is multiplied by an intensity, as it is in the equa- 
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tion of transfer, we get a term which can be neglected in comparison 
with the emission term. Consequently, 


dF jj ~ 
dt = een + By ji Mr dr . (56) 


Now, equation (48) shows that 


dF 


= > By Mii +. B, Ud, fo 4uiJ e 
dt — = 


k=j 


It is possible to solve these equations for 9. We suppose that 


dF 
Paar ~ 
Mi —, > je] dt ot 1 WU Sy . (57) 


To determine the coefficients, introduce this into the previous equa- 


tion: 
dF; am cle dF | 
=—_ -_ “N ( ° XN # So Pes \ r 
dt eg Pri ia dt O): Ur Dr oa Bj Mrdr 
k=j l=k | 
r c 
a dF; zs 
° + ei dt = Wj Ure , 
b=y 
~ dF Bee ~ - | . ‘is ~ » | 
as af dt } ~— J Br joi f - Udy ‘ Be; — Wy + Bij fe 
i oe k=j 
Therefore, 
co Wj; Bj jj; ’ 
1 
‘ —*, 
Oo = @j) — ‘\ Bre jWet , (1 > 7) 
k=j 
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From these follow the equations: 








l r 
~ PF aes = 
S Br j@ 41 Bri — S Bi iW 
ms ™ pie! BE > k=jtt 
Oo; = an = ———— oO; = ; 8 
i ae Bij’ " ies Bij ‘ ’ | a Bi; (5 ) 


Equation (57) can now be introduced into (56) and the integration 
performed. The result is 


Fi = — DY Dd BuiouF + (5? — 3)( + Yb 
k=j l=k k=j 


Using (55) and (58), we can write this in a slightly different form: 


F ij — _ | F; =e Sd aiF + wo; 3? (1 ran eet) ’ (59) 
~% l=j 


where we have substituted for 6;; its value p;,. Here we have utilized 
the vanishing of the fluxes on the inner surface of a planetary nebula. 

As an alternative procedure for these transitions we could have 
considered the specific intensity directly. For the equation of trans- 
fer we have 


addi; 


cos 6 
dt 


SY) Beiiuu( Je + 43:) . 


k=j 


Since the right-hand side is a function only of ¢, we can integrate 
this equation as it stands. Evidently the intensity will be propor- 
tional to sec 0, so that it will increase very rapidly with increasing 
values of the angle of observation. The actual infinite value of the 
intensity must not be taken seriously, since it is due to the approxi- 
mations involved—the neglect of absorption effects and the repre- 
sentation of the nebula as a plane-parallel sheet. Struve? has found 
that in certain diffuse nebulae the conditions necessary for the ex- 
istence of the effect are satisfied. His explanation is essentially the 
same as that given here. 


7 Ap. J., 85, 194, 1937- 
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This completes the formal development of the theory for a set of 
discrete states. The remaining steps are those necessary to general- 
ize the equations to the continuum. 


THE CONTINUUM 

Ambarzumian® has shown that the continuum can be treated ap- 
proximately as a discrete state, provided the weight is chosen ap- 
propriately. More recently Menzel and Pekeris’ have shown that 
each element of the continuum can be treated as a discrete state, 
when the weight of the element is taken according to the formula of 
Fowler. Now the number of atomic systems in velocity range (2, 
v + dv) for the free electrons is 


nif(v)dv, 


where m; is the number of ionized atoms in unit volume and f(z) rep- 
resents the distribution of velocities. The statistical weight of this 
element is 

8rmiv7dv 


nhs” 


where n, is the number of free electrons in unit volume. Hence, the 
number of systems per unit volume per unit weight is 


h3 NiMef (v) 


8arm3 v? 


NC v) = 


Without additional comment it is evident that the equation of 
transfer for the ith continuum can be put in the form (cf. [ro], [11], 
and [12]) 


dQ; (v) 5 3 7 ! 
1 = BONO) + (NO) — NJ uM), (12’) 


where the upper state is denoted by the frequency corresponding to 
the electron velocity, that is 
mv? 


2h ° 


8 Op. cit. 9M.N., 96, 77, 1935- 
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In (12’), 2:(v) represents the photon intensity in unit frequency in- 
tervals. It must be emphasized that the quantities B;(v) and yu;(v) 
are atomic functions and consequently do not vary from point to 
point in the nebula. 

Strictly speaking, (12’) gives only the partial variation of the in- 
tensity due to transitions between the states indicated. In general, 
the various continua overlap each other and the continuum formed 
by free-free transitions. Actually the effect of this overlapping is 
quite small for the first continuum and is purely formal for the 
others, since for 7 > 1 (and for free-free transitions) the second term 
on the right-hand side of (12’) can be neglected, as we shall see, 
removing all coupling between the different types of transitions. 
For the present, however, we shall regard (12’) as giving the net ex- 
change between the two states. 

Similarly, 


dQ(v, v’) 


1, = BO, YN!) + INGO’) — NO) ul, »DO, »') (12") 


gives the partial variation of the intensity due to free-free transi- 
tions, or the net exchange of systems between the states v and v’. 

The complete set of equations representing the condition of radia- 
tive equilibrium is as follows (cf. [14]): 


i—I co ) 
> (NB; + (N; = Nj) mii fii rr 3 LV Bi; + (NV; sii Ni)mii fii} | 
j=l j=i+i1 
=f 1XODBO) + IV) = Nano JO} de" = 0, | 

DV OBiO) + VO) — Nida) J.) thee 
1=% | 
| 


+ {ewe v) + IN(y) — NO’)Jur’, ») PO’, v) dv’ 


| 
: | 
= {NO DBE, ¥) + VO) — NOMnG, Jl, vd! = 0 


For dilute radiation the solution of these equations is the generaliza- 
tion of (36) or of (44). 
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If we define the quantities p(v); and p(y, v’) by the equations 
Bi (v) 


p> Daal +f B(v’, v)dv’ 


B(v’, v) 
\ B..(v) +f B(v’’, v)dv’’ 


q(v)i = p(v)i + > Q(v)c pri + { q(v, v’) p(v’);dv’ , 
k=i+1 ils (43’) 


qv, »’) = plv, v!) + if gv, »””) po", »!)de” 





pty): 


( (41’) 





p(y, v’) 


’ ‘ 


and q(v); and q(v, v’) by 


the generalization of (44) is 


Vi ~~ : 
Vy. Bi: = = Mri fri +S ae * quibaygiu +f qv’ )j Mr ( v ‘Vv dv’ 2 | 
7 pet ine 
' (44’) 


| Be) + L Be’, Hit | iia 


+f qv’, v) p(y’) Y.(v’ dv’ | 


Equations (44’) follow from the same line of reasoning that we em- 
ployed in deriving (44), and therefore depend on the assumption of 
the high dilution of the incident radiation. Furthermore, (44’) jus- 
tifies our deductions concerning the effect of overlapping. 
Hence, in the case of stratification in parallel planes, (14’) leads 

to the generalized flux integrals 

ar: 

j= 


I 


Fiji == Fi -{~ Fi(v)dv = = @, 
eel 


I j= mad 


> (18) 
5" S50 + [50 v)dv’ -(~ Fv, v’)dv’ = a(v) | 


j=l 
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It should be remembered that the continuous fluxes in (18’) have a 
meaning only if we are able to separate the component continua in 
the blended spectrum. Without any difficulty it can also be shown 
that (18’) leads to the generalized K-integral 








Sf S Kie 18 2 % Kile) — HOD) | | 
et Se ;' ae > (29! 
tf @ Kv, ey 4 Flv v’) x (v, v’) 
pnt As ae ———— . dy'dv=o0. 
u(y, v’) } 





Finally, we remark that the generalization of the r-integrals is 
trivial and need not be considered further. 

The equations (44’) being available, it is possible to proceed along 
the lines followed in the preceding section in solving the differential 
equations of the radiation field. It is evident that the formal ex- 
position will differ in no way from that followed in the discrete case, 
and we shall arrive at generalizations of (50), (51), (52), (53), and 
(54), in which the summations have been extended to integration 
over the continuous states. In particular, (51) and (52) become the 
pair of equations 


j es ) 
ae sey sage, ~ v’)dv | 
a;(v) = Sati dv)? | ) , By); mr ix i(v ’) BC V, ")dv’ ps ) (51 


yf ) 
a(v, v’) = _— | B(v’, v) -{ a(v, v’’)B(v’, v’’)dv”’ | » (52’) 
of which the second is an integral equation. In addition, B(y, v’) 
depends on the solution of the integral equation in (43’). The dis- 
cussion of these integral equations is beyond the scope of the present 
investigation, and we shall leave the problem in its present stage. 
The general solution of the continuous problem is under considera- 
tion and will be reported in a future communication. 


FOUR-STATE PROBLEM 
We close our discussion with a consideration of the four-state 
problem for a planetary nebula. From (43) we have 


das = Pass Jaz = Paz t Pass 5 Y32 = P32; 
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also 


By = Pus By = PasP3 5 By. = Px + PasP32 5 B83; = p; ; 
B32 = p32, Bu» =1 


From the definition of \;, we have 








4 = 3u3(1 a P41) ) N= mt — P31) ) AZ = 0. 
Furthermore, from (51) we have 
dy, = : PazsP3 a ey a ere 
(=) (0 = Ps) oo i Pax) 
M4 
and from (52) 
.— om ) a= Pubs — fae , a=-—I1. 
my I M; 
’ 3") (o~ a) 2 
M4 


If we associate the fourth state with the continuum, the quantity 
u;/m, Will be a large quantity of the order of 10. Taking this value, 
we have 
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where we have made use of the fact that \,/A, is large. Now, except 
for small values of X,/,, the optical thickness in the continuum, this 
becomes 


3 
F(t.) = 2 X ro~# PuPss_®s_ gp) 4 agore—met 
34 ) P32 3h; + 2h, ( sf ) 2745 


which is a small quantity. Finally, 
F(t) = Sr — eo") — F(t) — F(t) ~ Sa — e") — F,(4,) « 
Now we have 


= I J Pu ~ PasP31 
Oy = ’ On = ’ Gg = = 
Pees Pa = Pur p:2(1 = Pat) 


Hence, 


5. a 5 Daa Sa a aaa 

In Figures 1-4 are shown the curves for the various fluxes as 
functions of wl, = 7. Various combinations of the p;; have been 
used in computing these. Perhaps the most striking feature is the 
weakness of ¥,, which is roughly one ten-thousandth of the other 
intensities. This is certainly independent of the fact that we have 
considered only four states, and it will be true in a more general dis- 
cussion. Hence, we are led to the fact that the Lyman series, except- 
ing Ly a, is much weaker than the other hydrogen series. The fea- 
tures of the curves for %,, and &,, are, of course, identical with those 
of Chandrasekhar."° 


YERKES OBSERVATORY 
March 1938 


10 Op. cit. 











SODIUM IN THE UPPER ATMOSPHERE 
J. CABANNES, J. DUFAY, AND J. GAUZIT 


ABSTRACT 

Spectroscopic and interferometer observations show conclusively that the bright 
yellow radiation of the night sky coincides with the two D lines and that sodium atoms 
are always present in the upper atmosphere. At twilight the yellow radiation is in- 
creased, like the red line \ 6300 of [O 1]. The writers are examining the altitude of the 
luminescent layer and the origin of atmospheric sodium; they believe that the sodium 
comes from meteorites. 

I. THE YELLOW RADIATION OF THE NIGHT SKY AND 
ITS ATMOSPHERIC ORIGIN 

Among the lines or bands discovered by Slipher’ in the spectrum 
of the night sky, the bright yellow radiation near the sodium D lines 
calls for special attention. 

While the green line \ 5577 and the red doublet AA 6300-6364 
have for several years been definitely identified with the lines 
'S, > 'D, and 'D, > 3P.,, of the oxygen atom [O 1], the origin of the 
yellow radiation was still unknown a few months ago. 

The results of wave-length measurements taken by several ob- 
servers are summarized as follows: 


Angstrom 


: Units 
De NE i sow kein asses 1929 5892 
J. Defay?....... eT Te see 5892.5 
L. A. Sommerf. . b Peete eons gun) 5888 
De ee ee 1934 5888 
L. Vegard and E. Ténsberg||..... 10935 5885 
* Pub. A.S.P., 41, 263, 1920. § J. de Phys., 5, 601, 1934. 
t C.R., 194, 1897, 1932; J. de Phys., 4, 221, 1933. || Zs. f. Phys., 94, 413, 1935 


t Zs. f. Phys., 77; 374, 1935- 
The great discrepancies between these results have prevented earlier 
attempts to identify this line. New measurements were made on 
numerous plates obtained at Montpellier? in 1935, 1936, and 1937, 
using the Cojan spectrograph with large aperture, f/o.7. This 
spectrograph gives excellent reproductions of the yellow line at 

™ Pub. A.S.P., 41, 263, 1929. 2C.R., 206, 221, 1938. 
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twilight or at night, with half-hour exposures. Panchromatic “Ilford 
Hypersensitive”’ plates and Eastman F, Type I, plates were success- 
fully used. Near A 5900 the dispersion is about 500 A per millimeter. 
The width of the yellow radiation, on the plates, is 20 yu; errors aris- 
ing from the measurements of its center of gravity do not seem to 
exceed about + 1 um. 

The average value of the measurements taken on the ten best 
plates is 4 5894 + 1A. The yellow radiation of the night sky is 
therefore certainly quite close to the D lines. 

On the photograph reproduced here (Pl. VIa), part of the slit 
was illuminated by a sodium flame. The two D lines, which are not 
resolved by the spectrograph, have the same appearance as the sky 
radiation and are found in its continuation; however, they are per- 
haps displaced by 1 or 2 A toward the shorter wave lengths. 

In addition, the yellow radiation in our spectrograms has the 
same appearance as the D system and differs from the line \ 5577; 
the intensity of the yellow radiation decreases gradually on each 
side of the maximum, while in the case of the \ 5577 line the in- 
tensity decreases abruptly. The d 5894 line of the sky seems, there- 
fore, to be a complex radiation. 

A microphotometer tracing of the spectrum (the microphotometer 
enlarges seventy times) shows that the width of the line \ 5894, 
measured at mid-height, is greater than the widths of the lines \5577 
and \ 6300; the presence, on either side of \ 5896, of two lines a few 
angstroms away explains clearly this difference. The yellow radia- 
tion of the sky and the unresolved sodium D lines give similar 
diagrams. 

Observations by Dufay in 1932 were made by means of a fine 
wire set in the middle of a wide slit. The shadow of the wire was al- 
ways narrower and less clear in the yellow image than in the image 
of the green line \ 5577. This fact also seems to prove the complexity 
of the yellow radiation. 

Dufay examined, in 1932, the possibility of emission by inter- 
stellar sodium atoms; but the difference in intensity at the zenith 
and near the horizon and, in particular, the variations in intensity 
from night to night and during the same night gave evidence that 
the radiation was emitted in the upper atmosphere. Furthermore, 
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a long series of observations made by Garrigue at the Pic-du-Midi! 
put the atmospheric origin of the yellow radiation of the night sky 
beyond all doubt: the intensity of the radiation becomes two or 
three times greater from the zenith to the horizon. In addition, 
although the yellow radiation does not undergo regular variations 
during the night, as does the green line \ 5577 A, it seems to show 
fairly well-marked annual variations, the maximum occurring dur- 
ing the winter and the minimum in the summer. 

The hypothesis of an emission of sodium in the upper atmos phere, 
however strange it may appear to be, cannot be neglected. 


II. INTERFEROMETER ANALYSIS OF THE YELLOW 
RADIATION AT NIGHT 

The verification of this hypothesis requires the use of an ap- 
paratus able to resolve the two D lines. Since the spectrographs 
which we are using at the present moment have great light-gathering 
power but are of low dispersion and are not capable of separating the 
D lines, we have tried to analyze the structure of the radiation of 
the night sky with the help of an interferometer, which is easier to 
construct and is less costly than a prismatic apparatus. 

With a Fabry-Perot interferometer with air blades between half- 
silvered surfaces, the fringes to infinity resulting from the D, radia- 
tion insert themselves exactly between the fringes resulting from 
the D, radiation, when the thickness of the air blade is 0.15 mm. 
When the thickness is 0.30 mm, the two systems of fringes are, 
on the contrary, superposed. Working alternatively with these two 
thicknesses, we can, therefore, detect whether the yellow radiation 
of the night sky has the same structure as the D lines. 

The Fabry-Perot interferometers were composed of two plane- 
parallel glasses lightly silvered on one side and separated by three 
pieces of metallic wire. Invar wires 0.15 mm in thickness and 
platinum wires 0.30 mm in thickness were used. 

The first method.—In the first series of experiments the interferom- 
eter was adjusted in front of a small two-prism spectrograph the 
camera lens of which has a relative aperture f/1 and a focal length 
of 25 mm. A lens of 6 or 10 cm focal length projects the fringes on 


3 Ibid., 202, 1807, 1936; 205, 491, 1937. 
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the slit. With exposures of 20 hours or more, we have observed dis- 
tinctly, with the 0.15 mm interferometer, very crowded rings in the 
yellow radiation. The aspect of these fringes was exactly the same 
as that of the fringes given by the two D lines. At the same time, 
the fringes of the \ 5577 and Xd 6300 lines of [O 1] were intense. Agfa 
ISS panchromatic plates, hypersensitized before use in a diluted 
solution of ammonia, were used. The first positive result was ob- 
tained with a 27-hour exposure from February 28 to March 4, 1938.4 

The second method.—The intensity of the yellow radiation makes 
it possible to obtain interference fringes without the use of a spectro- 
graph, provided the \ 5577 and \ 6300 lines are eliminated. An 
orange filter (Schott glass OG2, 2 or 4 mm thick) can be used to 
absorb the green line. In order to be free from the red line, it is 
necessary to use plates whose sensitivity is weak at \ 6300. Cru- 
miére Panchro P plates fulfil this condition. 

We have photographed in this manner interference fringes at 
night with the Fabry-Perot etalon placed in front of a camera lens 
of 25 mm focal length and a relative aperture f/1. Using alternative- 
ly thicknesses of 0.30 and 0.15 mm, one or two systems of fringes 
appeared. The fringes are several millimeters in diameter and can be 
accurately measured. In every case their diameters differed by al- 
most o.o1 mm from the diameters of the sodium fringes photo- 
graphed on the same plate. The D, fringes are, as in the laboratory, 
about twice as intense as the D, fringes. Plate VIb shows the en- 
largement of a plate obtained from March 30 to April 2, 1938 (ex- 
posure 22 hours), with the 0.15 mm etalon. The copies of the 
fringes were cut along a diameter. We see that the fringes of the 
night-sky radiation are an exact continuation of those of the sodium 
flame. The sharpness of the outline, however, has been increased 
in the reproduction, and the former D, fringes are nearly indis- 
tinguishable on the sky photograph. One night’s exposure is suf- 
ficient to give distinct but faint fringes; clear photographs require 
two or three nights’ exposures. 

The preceding observations prove conclusively that the yellow radia- 
tion of the night sky is contained in the two D lines and that sodium 
aloms are always present in the upper atmos phere. 


4 [bid., 206, 870, 1938. 
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Ill. THE YELLOW RADIATION AT TWILIGHT AND 
ITS INTERFEROMETER ANALYSIS 


During the International Polar Year, 1932-1933, B. W. Currie 
and H. W. Edwards observed spectra of the polar aurora and of the 
night sky at Chesterfield, Keewatin, Canada.’ They found that the 
intensities of the yellow radiation and of the red line \ 6300 were in- 
creased during twilight.° However, their paper appeared only in 1936, 
and the strong excitation of the \ 6300 line had been already dis- 
covered and carefully studied by Cabannes and Garrigue at the 
Pic-du- Midi.’ 

The excitation of the yellow radiation was found again and 
studied by Bernard at Troms¢ in 1937." From his photometric 
measures it appears that the luminescence ends abruptly when the 
solar rays tangent to the earth reach an altitude of about 60 km. 
According to Bernard, there is a thin luminescent layer at this alti- 
tude, and the yellow radiation is excited by optical resonance. 

The interferometer analysis of the yellow radiation at twilight is 
much easier than at night, because of its stronger intensity. We have 
carried out the same experiments at twilight as at night. By our first 
method we photographed, with the spectrograph, fringes quite simi- 
lar to the sodium fringes, with thicknesses of 0.15 and 0.30 mm in 
the interferometer. Exposures of 1o-20 minutes are now sufficient. 
By the second method we regularly photographed interference 
fringes with exposures of 20-35 minutes. The exposures began in 
the evening when the sun was about 8° below the horizon. The 
camera was pointed toward the west, 10° or 20° high. In the morning 
the camera was pointed toward the east and the exposure was closed 


5 Terrestrial Magnetism, 41, 205, 1936. 

6In a recent paper (C.R., 206, 1137, 1938) R. Bernard attempts to show that the 
observations of B. W. Currie and H. W. Edwards did not concern the \ 5893 radiation. 
The wave length measured by Currie and Edwards was, indeed, \ 5940; but their 
measurements were not precise and they concluded that “this radiation is due to a 
luminescence distinct from that of the aurora and is probably the radiation of 5893 A 
observed by Vegard and T¢nsberg in the night sky.’”’ Therefore, Bernard’s argumenta- 
tion does not seem to us convincing. According to Currie and Edwards, the intensity 
of the yellow radiation is also increased during periods of moonlight. 


7 Ibid., 202, 1807, 1936; 203, 484, 1936. 
8 [bid., 206, 448, 1938. 
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when the sun reached the same altitude. The fringes were also 
photographed at the zenith, later in the morning and earlier in the 
evening, with a 15-minute exposure. 

Plate VIc and d shows enlargements of plates made at twilight 
with the interferometer 0.15 and o.30 mm thick. The identity in 
appearance with the sodium fringes is very noticeable; and, in the 
case of the 0.15-mm thickness, the intensity ratio of the two com- 
ponents seems to be the same in the sky as in the laboratory. 

At the same time, Bernard published quite similar results ob- 
tained during twilight, following the second method.’ 

In the conclusion of this article we shall deal only with the emis- 
sion of the yellow sodium radiation during the night. 


IV. RELATIVE INTENSITIES OF THE TWO COMPONENTS 
OF THE DOUBLET 
Our first measurements placed the center of gravity of the doublet 
at 5894 + 1 A. We concluded from this that in the night sky the 
two components D, and D, had the same intensity (as is the case in 
the spectra of comets"). The interferometer studies have not con- 
firmed this result: Plate VIb shows that component D, gives fainter 
rings than D,; the intensity ratio is near the theoretical ratio 1/2. 


V. ALTITUDE OF THE LUMINESCENT LAYER 


The observations made by Garrigue at the Pic-du-Midi in 
December, 1935, and January, 1936, give B,/B, as the ratio between 
the intensity of the yellow radiation at the horizon and at the zenith. 
The average of the intensities, which are very close to each other, 
as measured by Garrigue in the different azimuths at about 10° 
above the horizon, is 2.8 times greater than the intensity at the 
zenith. From this number we may derive the altitude of the lumi- 
nescent layer, supposedly homogeneous and relatively thin. When 
we observe one point in the sky, the apparent brightness is the sum 
of a term TBy, which gives the real brightness B, of the sky at that 
particular point, reduced by atmospheric absorption, and a term B’, 


9 Tbid., p. 928, 1938. 
© A. Adel, V. M. Slipher, and R. Ladenburg, Ap. J., 86, 345, 1937. 
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which is the result of the scattering by the lower atmosphere, for 
the whole length of the visual ray, of the light given out by the 
entire sky. The equation 


B= TB,+ B’ 


may be applied to the following four cases: the monochromatic 
luminescence of the sodium layer at 10° above the horizon and at 
the zenith, the scattering of the light of the stars by the interstellar 
particles, following the same directions. 

The observation of the continuous spectrum of the night sky gives 


Tat oo 

Bu B, , 

_ | wo. (1) 
1 z + B, 


An a priori calculation" shows fairly exactly the light diffused along 


a vertical ray: 


B! 

—- = 0.03 ; 

B, 3 ’ 
or a consequence of formula (1): 

By 

> = 0.105. 

B, ° 


This quantity is required to enable us to derive from the observa- 
tions of the yellow radiation the altitude of the luminescent layer. 
Let / be this altitude. The brightness of the layer at the zenith dis- 
tance z is the product of the brightness of the layer at the zenith and 


"This calculation has been carried out previously by J. Dufay, thése, p. 133, 


Paris, 1928. 
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of sec a, with sina = (R/R + h) sins (R = the radius of the earth). 
We have, then, in the case of the luminescent layer, 


en By 
Bu _ T sec aa 
\ B! 


——— Hy ots 





From this we conclude that sec a = 3.83. 
The altitude of the luminescent layer would then be about 130 km. 


VI. FREQUENCY OF THE TRANSITIONS *P > 7S PER SECOND 

The relation between the intensity of the yellow doublet and 
that of the green line of the aurora changes continually from one 
night to another and even during the same night. During the winter 
of 1935-1936 it varied from 0.85 to 1.3. But if we take as unity the 
average intensity of the green line, the average intensity of the yel- 
low doublet has the value 0.48. Taking the value 1.35 erg/sec/m? 
for the average energy radiation of the green line in the south of 
France at the beginning of the winter of 1935-1936, we get for the 
yellow radiation 0.65 erg/sec/m?’, this being the energy released by 
2.10'' transitions of the type 7P > ?S. 


VII. ON THE ORIGIN OF ATMOSPHERIC SODIUM 


At 130 km above the earth the temperature of the atmosphere 
does not increase beyond 200° C. We cannot believe that a thermal 
effect explains the emission of the D lines. The problem is, more- 
over, twofold: What is the origin of the atmospheric sodium, and 
why does it become luminous in the upper atmosphere? We can 
examine whether the sodium atoms have a terrestrial or a cosmic 
origin. The first alternative seems to be untenable: we do not under- 
stand how a fine dust of sodium salt can be dissociated in its ascent 
from the ground to 130 km. 

We tend to support the hypothesis of a cosmic origin. The earth 
receives each year from meteorites a deposit of 4 grams to the square 
kilometer. The proportion of sodium in weight is 7.2 X 10-3. We 
receive then, per second and per square meter, 2.5 X 107 atoms of 
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sodium. It may be supposed that the yellow radiation of the night 
sky is a luminescent phenomenon accompanying the fall of the 
meteorites; in this case each sodium atom would undergo eight 
thousand transitions before definitely dying out. 

The presence of sodium in the higher atmosphere, whatever may 
be its origin, leads us to look for other elements. In the stony 
meteorites, as in terrestrial rock, the most numerous atoms are” 
(after oxygen atoms) Si, Mg, Fe, S, Al, Ca, and Ni. We cannot ex- 
pect anything from Si, Mg, and S atoms; the Fe lines are too numer- 
ous and consequently too weak. Only calcium and aluminum must 
give lines with marked intensity. The resonance line of the neutral 
Ca atom has the wave length 4226.73 A; we have actually observed 
a line of \ 4226 in the spectrum of the night sky. It would be inter- 
esting to investigate whether the bright emission lines of aluminum 
do not coincide with certain lines of the night-sky spectrum. 

UNIVERSITE DE PARIS AND OBSERVATOIRE DE LYON 
May 1938 
27, and W. Noddack, quoted by V. M. Goldschmidt, Fortsch. d. Mineralogie, 19, 


153, 1935- 














ON THE INTENSITY DISTRIBUTION IN THE 
BANDS OF COMETARY SPECTRA 


P. SWINGS AND M. NICOLET 


ABSTRACT 

The exceptional intensity distribution of bands of unsymmetrical molecules like 
CN and CH ina comet spectra can be explained by taking into consideration the follow- 
ing factors: (a) the frequency of electronic absorption processes; (6) the distribution in 
the vibrational levels (ground state) which is a function of the nuclear temperature 7, 
or of the equivalent radiation temperature 7); (c) the distribution in the rotational 
levels corresponding to the nuclear temperature. These factors are different functions 
of the heliocentric distance r. 

1. Important progress has been made in recent years in our 
knowledge of the physical mechanisms which take place in the 
atmospheres of the comets; the most outstanding contributions in 
this field are those of K. Wurm.’ It seems useful to undertake a 
systematic rediscussion of the identifications of the emission bands 
observed in the spectra of comets; such a systematic work has been 
started in the department of astrophysics at Liége, Belgium.? 

Actually, very great care must be taken in such identification 
work, as the intensity distributions and the structures of cometary 
bands differ widely from what we observe in normal laboratory 
spectra. This point has been advocated by Wurm on several oc- 
casions. In the course of our systematic discussions we had the feel- 
ing that there were still several points to be cleared up in connection 
with the prediction or discussion of the intensity distributions, as 
well as with the physical mechanisms present. The aim of this paper 
is to bring some precision into certain of these points. We shall re- 
strict our considerations to the intensity distributions in rotational 
lines and vibrational bands of molecules which possess an electric 
dipole in their lowest electronic and vibrational state (such as CV 
or CH, but not C,); the fundamental difference between the be- 
haviors of the two types of molecules has been discussed by Wurm. 

' Cf. Zs. f. Ap., 15, 115, 1938; Handbuch d. Ap., 7, 305, 1936. 

2 The first contribution of our department in this field has appeared recently; see 
M. Nicolet, ‘‘Les Bandes de CH et la présence de l’hydrogéne dans les cométes,”’ 
Zs. f. Ap., 15, 154, 1938. Other papers by Nicolet are in preparation. 


173 











174 P. SWINGS AND M. NICOLET 


Owing to the very low pressures prevailing in cometary at- 
mospheres, the mean free path of the molecules is considerable. On 
the other hand, the excitation of the molecules has its primary 
origin in the solar radiation, the excitation by collisions being ex- 
cluded. Under these conditions, what is the mechanism which pre- 
vails in the distribution of the molecules in the rotational levels of 
the lowest electronic and vibrational state? 

2. Let r be the distance, in astronomical units, from the comet to 
the sun. By the application of Stefan’s law, and assuming To = 
5740 K, it is found immediately’ that a body absorbing the solar 
radiation at a distance r would take on a temperature4 7, ~ 3007 '””. 
But this factor, 7,, has no immediate meaning for the population in 
the discrete molecular levels in a gaseous atmosphere where collisions 
are practically absent. 

3. We may, furthermore, try to find for different wave-length 
regions some kind of equivalent temperature 7 in complete similar- 
ity with the case of interstellar matter;> this equivalent temperature 
T, must be such that the actual energy-intensity for wave-length \ 
is equal to that of equilibrium radiation at temperature 7,. This 
treatment seems at first logical in this problem, owing to the selec- 
tive absorption between the different rotational levels. Assuming 
that the sun radiates like a black body, and using the appropriate 
dilution factor, we find for the distribution in the excited levels, at 
the distance r = 1: for Selective absorption at \ ~ 4500 A, the 
equivalent temperature 7, ~ 2000° K; for selective absorption at 
A ~ In, the equivalent temperature T, ~ 800°; for selective ab- 
sorption at A~ 10 y, the equivalent temperature 7, ~ 70°; for 
selective absorption at }\~ 100u, the equivalent temperature 
T, ~ 5°. In any case, since the pure rotational lines of the mole- 
cules with which we are concerned lie in the region from too to 
1000 w or more, we should find by application of Boltzmann’s law 


3 See, for instance, A. Adel, V. M. Slipher, and R. Ladenburg, Ap. J., 86, 345, 1937. 
4For TO = 600¢°, T; would be 316r~'/?, which means a difference of 16° at r = 1. 


5 A. S. Edington, /nternal Constitution of the Stars, 2d ed., p. 371; for an application 
to the rotational levels of possible interstellar molecules, see P. Swings, M.N., 97, 212, 


1937. 
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that all the molecules are practically in their lowest rotational level. 
This does not correspond with the observations. 

4. Obviously, such a concentration in the lowest rotational level 
would be maintained only in the absence of perturbing processes. 
Collisions do not play any real role. But Wurm has called attention 
to an effect due to the electronic absorptions and emissions influenc- 
ing the distributions in the rotational levels. He has shown that, 
owing to the differences in intensity of the R and P transitions in 
emission and absorption, higher rotational quantum numbers may 
be reached, as long as the frequency of the absorption processes in 
the photographic region is higher than 1/7,, 7, being the lifetime of 
the rotational levels. 

This mechanism presumably takes place, but it is not easy to pre- 
dict its actual effect. On the other hand, the CV and CH molecules 
show a very different distribution; the CV molecules have their 
maximum population in the rotational levcls near A” = 8 (atr = 1), 
whereas the CH molecules are concentrated in the two or three low- 
est levels.© With Wurm’s mechanism, the observed difference brings 
complications and probably implies quite different 7, for CN and 
CH. We shall see later that, on the contrary, this difference appears 
conclusively when we assume another type of distributing mecha- 
nism. 

5. It is possible to reach a result quite similar to that of Wurm 
by an alternative theory which enables an easy calculation of the 
effect in the general case; in fact, there are several other reasons 
which support our hypothesis and which will be developed in the 
following paragraphs. 

The absorption of the solar radiation brings the cometary 
meteoric nucleus to a temperature 7, given by the application of 
Stefan’s law, following section 2. In the absence of any perturbing 
process the rotational populations of the lowest electronic level of 
any molecule will be practically the equilibrium distribution’ cor- 


6 Following the identification by Nicolet, loc. cit. 

7 As a first approximation we may neglect here the dilution factor; thus, the equiva- 
lent temperature due to the nucleus is the same for all wave lengths. On the other hand, 
it is obvious that the frequency of transitions between rotational levels in equilibrium 
at approximately 300° is much higher than at 7, ~ 5° or less. 
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responding to 7, = 300r*’?. The influence of the distance 7 is thus 
apparent. 

In certain cases Wurm’s mechanism may complicate somewhat 
the distribution in the lowest levels; but we may state that a rota- 
tional distribution corresponding to 7, = 30cr~'’’ seems sufficient for 
the interpretation of the observed features, as will be seen later.* 

6. We shall immediately show how simply our assumption ex- 
plains the different distributions of the CV and CH molecules. These 
molecules have very different moments of inertia: 


~ 


"(CN) 
I''(CH) 


14.10 ” gm cm ; 
2.10 ” gm cm’. 


Introducing these figures into the expression of the rotational energy 
in the Boltzmann factor, we find that they explain immediately the 
difference in distribution. 

7. The problem may be different in the case of the vibrational 
levels. Let us assume a frequency of vibration w’’ = 1000 cm™. 
The first vibration-rotational bands will lie between 10 w and 1 yp; 
in this interval the equivalent temperature 7,, following section 3, 
ranges from 70° to 800°. Thus, for the lowest levels the ‘‘nuclear”’ 
temperature 7, will have the predominant effect, but not for the 
higher-lying vibrational states. 

8. Summarizing the preceding considerations, we may say that, 
excluding intermittent of exceptional phenomena, the intensity of a 
cometary band depends upon the following factors: (a) the fre- 
quency of electronic absorption processes, which depends upon the 
spectral region and the structure of the absorption bands; (b) the 
distribution on the vibrational levels of the electronic ground state, 
which is a function either of the nuclear temperature 7, (for the 
low quantum numbers 2’) or of the equivalent temperatures 7, (for 
high v’’); (c) the distribution on the rotational levels corresponding 
to the nuclear temperature 7, = 300r-'’*. These factors are different 
functions of the distance 7. 

g. A verification of our hypothesis may be found in the good 
agreement obtained when we apply it to the identifications of CV 

* An elementary calculation shows that the frequency of pure rotational transitions 


of molecules in equilibrium at 7; = 300° is at least as high as the frequency of electronic 
absorptions near \ 4500 for an equilibrium at 7, ~ 2000”. 


BANDS OF COMETARY SPECTRA 177 


or CH/ lines or groups of lines.’ We shall not give here all the details; 
these will be published by one of us (M. Nicolet). 

For any identification the resolving-power of the spectrograph 
used’? has to be considered first; the theoretical intensity of the 
rotational line must be examined next. Let A be the width in 
angstrom units of the monochromatic image of the slit on the plate. 
A convenient method for a first approach will be to represent each 
rotational line by a rectangle of width A and of height J, J being 
the theoretical intensity of the line. J is obtained immediately by 
classical formulae introducing the so-called intensity factor 7 and 
the Boltzmann expression. If A is smaller than the wave-length 
separation between two lines, these will appear separated on the 
spectrogram; but in the opposite case, the intensities will be dis- 
torted and may be represented by the sum of the rectangles which 
are superposed. Thus, we may find a blackening on the spectro- 
gram, either (a) where the lines are very close together and not too 
weak (in the neighborhood of certain band heads) or (6) for the rota- 
tional quantum numbers giving maximum values of J. It may thus 
happen that a single branch will exhibit two marks on the plate: one 
near the head and one for the lines corresponding to the maximum 
values of 7. For instance, the P branch of the (0, 0) transition of 
CN, which has a head, may show two maxima; inversely the R 
branch, which has no head, will exhibit only one maximum. In fact, 
this R branch is blended with the P branch of the (1, 1) transition; 
thus, when the emissions starting from v = 1 are strong enough, 
this blending effect will complicate the aspect of the R(o, o) branch, 
although it is possible to interpret all the observed features." 

10. The main results concerning C'V may be summarized as fol- 
lows: 

Let us take A = 2 A, which is approximately the case for the 
spectrograms of the Comet Brooks taken by W. H. Wright.” As- 

9 We may safely say that many identifications, even of molecules like CN, are in- 
complete or erroneous! 

© We shall be concerned here exclusively with slit spectrograms. 

‘1 We want to emphasize the fact that the appearance of the spectrograms of the 
same comet obtained at the same moment with instruments of different resolving- 
power may be very different. See W. H. Wright, Lick Obs. Bull., 7, 8, 1912 (Comet 
Brooks). 

2 Tbid., n. 11. Pl. No. 337: A = 1.9 A; Pl. No. 338: A = 2.4 A. 
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suming for T,, respectively, 280° (r ~ 1, 2) and 420° (r ~ 0, 5), we 
have drawn diagrams of the intensity distributions in the (0, o) 
bands in two cases: case A: without considering the effect of the 
(1, 1) transition; case B: with consideration of the effect of the 
P(1, 1) branch on R (0, 0). These diagrams are shown in Figures 
1, 2, and 3. The following results appear at once: 

a) The band head's near \ 3881 has a rather sharp edge on the 
long wave-length side. Wurm has tentatively explained this phe- 
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lic. 1.—P and R branches of the (0, 0) transition in CN. 7; = 280° K. Without 
consideration of the P(1, 1) band. 


nomenon" by considering.the intensity of the solar radiation and the 
density of the solar absorption lines in the exciting wave-length 
region. The effect mentioned here gives an alternative explanation ; 
it is particularly conspicuous at T, = 420° (see Fig. 2). 

b) Beside the head at A 3881, there is another maximum near 
AA 3878-3879."5 The head will be prominent at small distances 7 
(T, ~ 420°); the opposite will happen at large distances 
(T, K 300°). 

c) The band near \ 3870 is due to R(o, o) in case A and to the 
blending of R(o, o) and P(1, 1) in case B (Fig. 3). Case A corre- 

13 The laboratory head is at » 3883.5. 

ZS: AD: 8, 10; 1932: 


‘STf A < 2, this maximum may be separated from the band head; if A > 2, the 
separation will disappear and the microphotometer structure will be different, the 


principal feature being either the head or the maximum. 
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sponds to great values of the distance r (low population of the vibra- 
tional level v’’ = 1); Figure 1 gives us, then, a diffuse band with 
maximum at AA 3869-3870. Case B will happen for small distances 
r; the maximum of case A at d 3870 will then be distorted by another 
maximum of the P (1, 1) branch at \ 3868. 

d) Thus, at a very small distance r we shall find simultaneously in- 
tensity maxima at A 3881(head of P[o, o]) and \ 3868 (maximum of 
































Fic. 2.—P branch of the (0, 0) transition in CN. Continuous curve: 7; = 420 
(r ~ 0.5); dotted curve: 7; = 280° (r © 1.2). 

Fic. 3.—R branch of the (0, 0) transition in CN. Continuous curve: 7; = 420 
without blending effect of P(1, 1); dotted curve: with consideration of P(1, 1), as- 
suming for convenience /(1, 1) = 3/(0, 0) (J = intensity). 


P(t, 1]); at a very large distance, intensity maxima will appear at 
3879 (maximum of Plo, o]) and AA 3869-3870 (maximum of 
R\o, o}). 

e) Similar results are obtained for the bands near \ 4200, for 
which our assumption brings clearly an intensity-distribution quite 
identical to the (0, o) band. The band head near \ 4215 is especially 
conspicuous for T > 420°, and the duplicity of the band at d 4196 
is due to the superposition of R(o, 1) and P(1, 2). 
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f) An examination of the wave lengths obtained with slit spectro- 
grams" discloses a satisfactory agreement with the preceding con- 
siderations; of course, we have to consider carefully the effects of the 
distance r. Similarly, there is no difficulty in explaining certain dif- 
ferences between the spectra of different comets. 

g) The subdivision of certain CN bands does not require special 
excitation processes of violent character; it is just a consequence of 
the preceding considerations. 

h) Similarly, there is no need for two types of excitation, as sug- 
gested by A. Adel.'? For Comet Zlatinsky (r ~ 0.7 a.u.) the band 
at \ 3881.1 corresponds, as has been stated by A. Adel and Wurm, 
to the wave lengths of the rotational lines P (13) or P (14) of the 
(o, 0) transition; but it is important to notice that actually the 
theoretically strongest line is P (8) or P (9); the observed effect is 
the result of the superposition of lines, owing to the value of A. 
This appears conclusively in the diagrams, as well as in the behavior 
with varying distance r. 

i) It is obvious that a Boltzmann distribution at temperature 7, 
is only a first, and possibly rather crude, approximation. The popu- 
lation in the lowest rotational level may be higher than in the case 
of a true Boltzmann distribution. If so, the lines observed near 
AA 3864, 3873, and 4190 could be, respectively, the origins of the 
(o, o), (1, 1), and (1, 2) transitions. 

j) A careful investigation of the red system of CNV in comets is 
still lacking, the actual data by F. Baldet,"® N. T. Bobrovnikoff,’® 
Dufay, Bloch, and Ellsworth,” J. Gauzit,?* and R. Minkowski” be- 

© Comet Daniel: Campbell, 4p. J., 28, 229, 1908; Comet Brooks: Wright, Lick 
Obs. Bull., loc. cit.; Comet Beliavsky: Bobrovnikoff, Pub. A.S.P., 43, 61, 1931; Comet 
Zlatinsky: Slipher, Lowell Obs. Bull., 2, 67, 1914; Comet Delavan: Curtiss and Mc- 
Laughlin, Pub. U. Michigan, 3, 203, 1923; Comet Mellish: Slipher, of. cit., p. 151. 
1916; Comet Finsler: Minkowski, Pub. A.S.P., 49, 270, 1937; Comet Halley: Bob 
rovnikoff, Pub. Lick Obs., 17, 443, 1931. 

7A. Adel, Pub. A.S.P., 49, 254, 1937. 

18 F. Baldet, thése, Paris (Annales de Meudon, 7, 1, 1926). 

19 Bobrovnikoff, op. cil., n. 16. 

7 Dufay, Bloch, et Ellsworth, C.R., 204, 663, 1937; J. de Phys., 8, 17, 1937. 

21 J. Gauzit, C.R., 206, 160, 1938. 


22 Minkowski, op. cit., n. 16. 
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ing only preliminary in character. Owing to the different electronic 
transition (*II > 72 instead of 7X > *Z), the intensity distributions 
among vibrational bands, as well as the rotational lines, must be 
quite different. 

11. The identification of CH by M. Nicolet” was based on the same 
principle. Adopting A = 2A and 7, = 300°, we obtain Figure 4. 
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lc. 4.— Intensity distribution in CH between \ 4290 and \ 4315 A. 7; = 300° 


This explains at once the following features: (a) \ 4313 is a rather 
narrow and strong emission; (b) \ 4297 and ) 4291.5 are blends of 
two rotational lines; the intensity of \ 4297 is similar to that of the 
single line \ 4304; (c) when A > 2 A, a single band is observed from 
4 4296 to A 4301 (blend of three lines). 

A discussion of the other electronic transition B *> > X *ILof CH 
will be published later by one of us (M. Nicolet). 

UNIVERSITY OF LIEGE 
DEPARTMENT OF ASTROPHYSICS 


April 1938 


23 Nicolet, loc. cit. 





THE ANALYSIS OF THE INFRARED LIMIT OF 
ATMOSPHERIC TRANSMISSION 


ARTHUR ADEL AND C. O. LAMPLAND 


ABSTRACT 

The correspondence between the related set of carbon dioxide bands at 15 uw and the 
infrared limit of atmospheric transmission is discussed in detail. A sample curve of 
atmospheric transmission in this spectral region is given by way of illustration of the 
detail actually present in what had formerly been considered a featureless, smooth 
curve, 

Description of the rock-salt prismatic solar spectrum in the neigh- 
borhood of the infrared cutoff—The appearance of the long wave- 
length end of the telluric spectrum may be briefly reviewed as fol- 
lows. It is characterized principally by the powerful absorption of 
ozone at 9.6 w, the cutoff in the vicinity of 14 uw, and the two small 
absorption bands gq, and g, between.’ It is the purpose of this paper 
to specify more precisely the character of the infrared limit of at- 
mospheric transmission. 

The infrared limit of transmission.—Figure 1 is a reproduction 
from the work of Martin and Barker of the related set of carbon 
dioxide absorption bands at 15 u.? The atmospheric opacity in ques- 
tion is produced by carbon dioxide by virtue of the existence of this 
set of bands.’ A high resolving-power curve of a portion of this re- 
gion, also from the work of Martin and Barker, is shown as Part a 
of Figure 2. Part b of Figure 2 is the corresponding section of the 
rock-salt prismatic solar spectrum . Part b was recorded at the Lowell 
Observatory on May 9g, 1938, during the interval 10:35-11:00 A.M., 
M.S.T. The sky was clear. The connecting lines in the figure indi- 
cate the corresponding regions of absorption. 

General remarks concerning the absorption bands.—The identifica- 
tion of the bands is well known and is given in connection with 
Figure tf. 

' A. Adel and C. O. Lampland, 4p. J., 87, 198, 1938. 

2 Phys. Rev., 41, 291, 1932. 

3A rather weak absorption band of ozone with its center at 14.2 4 contributes 
slightly to the cutoff. 
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Fic. 1.—The absorption of infrared radiation by carbon dioxide at 15 u 
IDENTIFICATION OF THE SET OF CARBON 
DIOXIDE BANDS AT 154 
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Fic. 2.—The infrared limit of atmospheric transmission 
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It can be shown that the matrix element of the electric moment 
for each of these transitions is of the order of 8, where Be'’ is the com- 
ponent of the electric moment perpendicular to the symmetry axis 


TABLE 1 


Approximate 
Wave-Length Interval Percentage 
Transmission 
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Fic. 3.—Atmospheric transmission at the cutoff 


of the molecule. The intensities of the absorption bands will there- 
fore be in the same relation as the molecular populations of the sev- 
eral initial levels. This fact explains the near-symmetrical nature of 
the laboratory-curves shown in Figures 1 and 2. 
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Concluding remarks._-The apparent percentage transmission of 
our atmosphere in this region of the spectrum can be obtained by 
extending the indicated envelope of the sun’s energy-curve. The re- 
sults are summarized approximately in Table 1 and in Figure 3 for 
the atmospheric curve shown in Part } of Figure 2. This curve of 
transmission will vary, of course, with the zenith distance of the 
sun. It may perhaps vary with the season and with the location of 
the observing station. 

LOWELL OBSERVATORY 


FLAGSTAFF, ARIZONA 
June 1938 








FURTHER DETAIL IN THE ROCK-SALT PRISMATIC 
SOLAR SPECTRUM 


ARTHUR ADEL 


ABSTRACT 

New features in the infrared telluric spectrum, made apparent with increased re- 
solving-power, are briefly discussed. 

Envelope structure of the ozone band at 9.6 u.—-With increased 
resolving-power it has been possible to disclose further detail in the 
rock-salt prismatic solar spectrum." 

We note in Part a of Figure 1, for example, that the envelope 
structure of the strong ozone band at 9.6 wu is clearly that of a doub- 
let. The character of the envelope of this absorption band is of 
considerable importance in deducing the structure of the ozone 
molecule and in analyzing the latter’s infrared spectrum.* The 
present prismatic measures, made at the Lowell Observatory, 
corroborate the grating measures made at the University of Michi- 
gan in 1935, in which this envelope form was first revealed.? The 
grating observations, which placed the band center at 9.57 uw, are 
shown in Part 6 of Figure r. 

Detail in the spectral interval from 7.2 to 8.5 w.-Of greater im- 
portance is the disclosure of additional absorption bands in the 
vicinity of the diverging wall of the great water band v,. The new 
curve is shown in Figure 2. Definite identifications are postponed 
for the future. As a preliminary, however, it may be instructive to 
point out the possibility that these new bands may indicate the 
existence of additional oxides of nitrogen in the atmosphere. For 
example, the molecule of V.O possesses two doublet bands with 
centers at 7.77 and 8.57 uw, the former being more than twice as in- 
tense as the latter.s It may well be that the R branch of the V,O 


' For a description of the telluric spectrum please see Ap. J., 87, 198, 1938. 
2 A. Adel, V. M. Slipher, and O. Fouts, Phys. Rev., 49, 288, 1936. 
3E. K. Plyler and E. F. Barker, Phys. Rev., 38, 1827, 1931. 
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FIG. 1. 


Envelope structure of the ozone band at 9.6 
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band at 7.77 u and the V.O, band at 7.6 uw combine to form the 
telluric absorption with its center at 7.63 y.' 
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Fic. 2.—Detail in the spectral interval from 7.2 to 8.5 u 


Detail d (yu) Detail d (us) 
wae , 44 fs 8.10 
b.. 7.55 g 8.17 
c 7.93 &r 8.30 
d 7.97 h 8.48 
ees 7.85 i 8.58 
: fe 8.05 j 8.72 


The wave-lengths are approximate. It should be remarked also that the deflec- 
tions at the base of the great water band are probably spurious. The detail enumerated 
in this table and shown in Fig. 2 has been observed regularly since its discovery. 


LOWELL OBSERVATORY 
FLAGSTAFF, ARIZONA 
June 1938 
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NONRADIAL OSCILLATIONS OF STARS 
PART I 
C. L. PEKERIS 


ABSTRACT 

In the first part of this paper a derivation is given of the Eulerian equations which 
govern the small and adiabatic oscillations of a general spherical harmonic type of non- 
rotating stars. These equations are then solved for a model of constant density. It is 
shown that with the exception of purely radial oscillations, which are unstable only when 
<4, there exist nonradial modes of oscillation, which are unstable for all values of y. 
The conclusion is drawn that in the pulsation theory of stellar variability it is necessary 
to consider nonradial disturbances in order to ascertain the possibility of the existence 
of any assumed model. 

I. INTRODUCTION 

In current theoretical discussions of stellar pulsation it is custom- 
ary to leave out of consideration nonradial modes of oscillation, with 
the justification that these modes are likely to be associated with a 
greater amount of viscous damping than purely radial oscillations. 
It may be argued, however, that in the lowest modes of nonradial 
oscillation the viscous damping is not likely to exceed by a great 
deal the damping in radial oscillations. Furthermore, we must con- 
sider the circumstance that, whatever external disturbances may 
exist, they will always be of a nonradial character. More serious, 
however, is the fact, pointed out by Rosseland,' that a limitation to 
purely radial oscillations fails to reveal the possible instability of the 
assumed model for a more general type of disturbance. If there is 
any mode of oscillation for which the assumed model is unstable, 
then the possibility of the existence of a star built on that model is 
nil; and therefore the periods of purely radial oscillations of the mod- 
el can have no application to actual stars. An investigation of non- 
radial oscillations should thus lead to an elimination of a large class 
of stellar models and would further add significance to the computed 
periods of radial oscillation of those models which prove to be stable. 
Among other problems which can only be solved by considering 
nonradial oscillations is, for example, the question of the degree of 


*S. Rosseland, Oslo Obs. Pub., 1, No. 2, 2, 1932. 
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instability of a star consisting of an adiabatically unstable core and 
a stable shell. 

In the following we derive the Eulerian equations which govern 
the small adiabatic oscillations of nonrotating stars. This problem 
was treated by Rosseland? in connection with his investigation of the 
effect of nonadiabatic processes on the stability of stars, but he made 
a special assumption concerning the functional form of the variable 
part of gravity. The variation with latitude and longitude of any 
mode of oscillation is given by a factor of a surface spherical harmon- 
ic, and we shall henceforth distinguish between the various modes by 
the order of this spherical harmonic. It is shown that the oscilla- 
tions are governed by an ordinary differential equation of the fourth 
order. In the special case of a model of constant density this equa- 
tion reduces to a system of two simultaneous equations of the second 
order, the solution of which is given. This solution yields, first, the 
known frequencies of purely radial oscillations, which are stable pro- 
vided y>4. The nonradial oscillations are, for any given spherical 
harmonic, either of a stable or an unstable type. The existence of 
these nonradial stable modes for a model of constant density is at 
first surprising, but an examination of one pair of modes for the 
spherical harmonic of the second order shows that in the stable case 
the destabilizing forces act in opposite directions in the central and 
outer parts of the star, thus neutralizing each other. In the unstable 
mode the motion is such that the destabilizing forces in the two parts 
of the star reinforce each other. It must be emphasized, however, 
that the existence of stable modes of oscillation is of no consequence 
if there also are other modes which are unstable, because under any 
type of disturbance, except of one class of extremely special nature, 
the unstable modes will grow exponentially with time relative to the 
stable modes and will eventually cause a transition from the origin- 
ally assumed configuration to another one. 


2. THE THEORY OF SMALL ADIABATIC OSCILLATIONS OF 
FLUID GLOBES 
The limitation to small oscillations implies that the amplitudes of 
the velocities and of their derivatives are assumed to be so small that 


abd 5. 22: 
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their squares and products can be neglected in comparison with their 


first powers. We also neglect viscous forces and assume that the 
motion of any individual particle occurs adiabatically. In spherical 
co-ordinates let w, u, and v designate respectively the radial, latitu- 
dinal, and longitudinal components of velocity. Further, let the sub- 
script “‘o” refer to the equilibrium value of any quantity, and let a 
symbol without a subscript stand for the perturbed part of that 
quantity. Then with 





the equation of continuity, 
es — ee (2) 
reduces under the foregoing assumptions to 
a. (3) 


ot or 


For an ordinary gas the condition of adiabatic motion is expressed 
by 


1 Dhotp)_ 7 Dio +) 


(+p) Dt (m+) Dt’ = 
where ¥ is the ratio of specific heat at constant pressure and constant 
volume. Under stellar conditions y is a variable quantity, but the 
limits of its possible variation are relatively small and we shall as- 
sume for simplicity that it is constant. It will also be convenient to 
introduce the quantity c, having the dimension of a velocity, defined 
by 
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In an ordinary gas, c is the Laplacian velocity of sound. Thus, (4) 
becomes 


¢) @ 0 d d 0 4q 
op sis Op = (2 (= +4. a, = — p,c*X. (6) 


© Or at | or 


ot Or 


Remembering that 


O Po 


ar — — £0 Po 5) 


we obtain from (6) 
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The gravitational potential will consist of an equilibrium part, y, and 
a perturbed part, y, which satisfy, respectively, the equations 


2 o_- G 094 ) ) 
sis 4m | (8) 
Vy —4nGp .)} 


Again, on neglecting quantities of the second order, the equations of 


motion become 


Ou °* 10 

Poa = Fo (p — po) , (9) 
Ov I 0 

ia at arg r sin 0 ap P — Po) » cael 
Ow _ Op ad 

po = ap + Po > — BoP - (11) 


We now assume that the dependence of each of the perturbed quan- 
tities upon time is given by a factor e'”, where the real part of o 
stands for (27/7), T being the period. Any possible mode of oscilla- 
tion will be characterized by an appropriate value of ¢, which will 
in general be complex; and the mode will be stable, unstable, or neu- 
tral if the imaginary part of o is, respectivey, positive, negative, or 
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zero. By differentiating (11) with respect to ¢ and (7) with respect 
to r, we obtain, with the aid of (3) 


ow — gw — gwtiop+cX+g(1 —y)X=0. (12) 


Now substitute for « and v in (1) from (9) and (10), obtaining 





wl _ (rw) re. : Re (s sin 6 =) + = ss = (é = ) . (13) 


or \sin 000 sin? 6 Og} \ po 


If the dependence of the perturbed quantities on @ and ¢ be given 
by a spherical surface harmonic of order 7, then the operator in { } 


is equal to —n(z + 1). Thus, with the aid of (7), (13) becomes 


ee aa = 


r? or or 





(gow — c?X) + Y, 





or 


rw + wi ar — golem + 1) 4 + n(n + oh 


(14) 
=x |r = n(n oS]. | 
o7| "| 
Finally, on differentiating (8) with respect to /, we obtain 
le Op ; abe 
ioVy = —4nG 5° = 4x (ws i x) (15) 


We now proceed to eliminate y from the system of equations (12), 
(14), and (15). First we compute V’*y from (12) and (14) and sub- 
stitute in (15), obtaining 


ve 4 26° | 
eX +X e-(y— e+ | 
; 2) | 
tei ”)(é + i — n(n + i > > (16) 


= go + wv ( 24 + = =) +l — set 1] = — 
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On differentiating (14), one obtains,’ with the aid of (12), 


Xr +X alg — 1g — le 2) S ends 
17 


= rwv+ arw + wl2 — n(n + 1)]. 


The elimination of w from the last two equations is rather cumber- 
some, and we shall give only the final result. Let 


A=090 g_t 1G = 2X ‘ iD ms y 2° 
A a(é =), g( ¢ +X (« ye + : ) 
+x [0° + (2 —yet in n(n + i 


n 





+ nt ye —g+8)], 


o*r? 


H=X+%[4-2 wt ne- st) 
r o*r 
r 2 n ee y . 
+x[2- et 1) (é ‘+70 ; 


then 


“ - (6 A x ee ( —-n-—-w 6A 2A? | 
ide, ag re Bae a 


r ia 





(18) 
— AH =o.| 


The boundary conditions are, in the case of a star of finite radius, 
first, that X shall be regular everywhere. At the boundary (r = R), 
D(p. + p)/Dt must vanish. But by (6) this is equal to —p,c?X or 
to —yp .X; and, since fp, vanishes there, the foregoing condition 
will be satisfied if X does not become infinite. We have an additional 
condition arising from the fitting of the perturbed gravitational poten- 
tial at the boundary. Let the dependence on r of the gravitational po- 
tential inside and outside the star be denoted, respectively, by y and 
y’; then, corresponding to any spherical harmonic n, the latter will be 


3 From now on we shall omit the subscript in go. 
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of the form (C/r"+"). The continuity of the gravitational potential 
at R requires that 
C 
Roti = Re (19) 
If, now, 7 denotes the amplitude of the deviation of the boundary 
from its equilibrium value Rk, we have the additional condition that 











Op dp’\ _ wt 
es: ade ) = 4mG(pon)r ; 
or 
Oy Oy’) _ —— 
a e re ia = 4rG(pow)p . (20) 
But 
oy’ 2 c . - VR 
(~) = —(n+ 1) on™ (n+ 1)p : 
so that (20) becomes 
. |d f 
io | + (n+ 1) we 47G(po.w)p . (21) 
or R R 


The foregoing system suffices to determine uniquely the state of 
motion of a star when perturbed slightly from its equilibrium con- 
figuration, as will be illustrated in one simple case in the next sec- 
tion. Before doing so, we shall verify that our equations agree with 
the Lagrangian equation of the second order, which is known to 
control purely radial oscillations. When m = 0, (14) reduces to 

2w 


X= w+ — (22) 


r 
and (15) integrates to 


lop = 4rGpw, (23) 
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which, when substituted in (12), yields the equation 


ba 7— —j= > 
o - + (2 53 oO. (24) 





= . | 20? 
Dw) = Cw + wh i? no v6 | +w 
By making the substitution w = rf(7) in (24), we obtain an equation 
for f: 

ey 7 (40 : io P 

of +7 ( - 18) +40 + (4-39 8] =o, (25) 
which is identical with the one in Eddington’s theory of radial oscil- 
lations. Equation (17) reduces to the result of multiplying (22) by 
r and differentiating once, while equation (16) reduces to 


d 2 

( agar Dw=o. 
dr r 

3. NONRADIAL OSCILLATIONS OF A MODEL OF 
CONSTANT DENSITY 


Let A = (47/3)Gp., 8 = (07/A); then 


bo |~2 


tf, = Ar, bo = 4Ap.(R? — 1’), c= — A(R? —-?r), 


and the right-hand sides of (16) and (17) become, respectively, 


\ 


rt fo tle =~ ne +a}. 
\ r w) 
Thus, on eliminating the expression in { } between (16) and (17), 
we obtain 

i. ae 2 8 2 
(R? — r?)X + — (2R? — 6r?) + X a —6—(n+1) Ae 
r ae 7B (26) 


-@+)2@-9)| = 0. 
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It will be convenient to introduce the nondimensional co-ordinate 
y = (r/R), whereby (20) transforms into 


— y)X += (2 — 6y’) 
= , (27) 
+ X e+ —6-@+)-@+)26—9)| =o: 


Having solved (27), ¥ can be determined from (15), which reduces 
to 
icpV*y = 3AX, (28) 


and then w can be determined by eliminating w between (12) and 
(14). Now, the solutions of X behave at the origin like r" or r-"*. 
Hence, we write X = y"F(y), and substitute it in (27), obtaining 


(1 - y)F be * [an + 2— y(2n+ 6)] + BF =0, (29) 


where B= — 6 —4n+ (8/y) + (28/y) — 2(n + 1) (n/yB). Let 
F= sv Cry"; then it follows from (29) that 
emma 


0 





; _ {| Ma+ s+ 20) - 28) 
“eee SE) + 2a + 3 + 28) | (30) 


When k becomes large, (¢x,2/c.) approaches unity, and it can be 
shown that F will not converge at y = 1. Hence, in order to make 
X(1) finite, we must have 


B= k(k+5+ 2n), k=O, 2, 4.2%; (31) 


and F will be a polynomial of degree k. Equation (31) may be writ- 
ten in the form 


B- (n+ 1) 5 = —4t FRR + 5 + 2m) +6 + an] = 2D, 
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or 
B=D+[(P+n(n4+1)}. (32) 


It is clear that, no matter what the values of k and y may be, there | 
is always one value of 8 which is negative, implying an unstable | 
mode of oscillation. When » = 0, the relevant root is | 


B= 2D =3y-44+7 kk+5), (33) 


as can be ascertained directly from (27). Equation (33) shows that | 
in the mode without nodes, for which k = o, B becomes negative 
when y < 4, a result which was first derived by Ritter. 

To obtain more insight into the nature of the modes which cor- 
respond respectively to the two possible values of 8 when n # o, we 
shall study the case n = 2, k = o. We have 


2D= 77-4, i= Cr, 


rr 
Hence, 
iop = ACri + Er’, (34) 
where C and £ are arbitrary constants. It follows from (12) and (14) 
that 
w ee — rg+ org — n(n+ 1) A = —crX 
o 


+X ls — n(n + 1)¢? a — iorp — in(n + 1)¥ ge 
0 0 
where the dots stand for differentiation with respect to r. In our 


case, (35) reduces to 


JB +1 — (w+ 1) lw = (ry — sw = — TR rik | 
| + (36) 


+X [a —o(n + nr) 7] i Fb - inet 5. | 
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The constant E in (34) can now be evaluated by substituting the 
value of we from (36) in (21), the result being 


E = — AACR’. 
Hence, 
Y= *ACr(r? — R’) (37) 
and 
aa 3 3 | 
w= CRIT (243 —iI- =). 8 
7B? +B) ~*~ 8] oe 


The expressions in { } for y = (5/3) and (4/3) are given in Table 1, 
where the superscripts (*) and (~) refer, respectively, to the positive 











TABLE 1 
Y pt hig B- ‘em 
C1 ee eee te 8.38 2.36(r?/R?2)—1. 36 —o.716 —2.109(r?/R?)+3.19 
NR irae, Ge eY ING 6.29 2.48(r?/R2)—1.48 —0.954 —1.14(r?/R2)+2.14 





and negative roots of (32). It is seen that in the stable solutions 
(positive 8) w changes sign in the interior of the star, while in the 
unstable solutions it does not. The stability of the former solutions 
is thus to be explained by the neutralization of the opposing de- 
stabilizing forces in the central and outer portions of the star. Sim- 
ilar results are found in the case m = 1, k = o. These modes are dis- 
tinguished from all the others by the fact that w does not vanish at 
the origin, which simply expresses the fact that a deformation of the 
form of a spherical harmonic of the first order is equivalent to a shift- 
ing of the center of mass. 

Further applications of the foregoing theory to other models will 
be made in a forthcoming paper. 


I wish to express my appreciation of the interest shown by Pro- 
fessor S. Chandrasekhar in this and in other problems of stellar 
hydrodynamics, and to thank Professor O. Struve for the oppor- 
tunity to spend some time at the Yerkes Observatory. 

YERKES OBSERVATORY AND MASSACHUSETTS 
INSTITUTE OF TECHNOLOGY 
June 15, 1938 
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ORIGIN OF THE INFRARED TELLURIC 
ABSORPTION BAND 4; 
ABSTRACT 


It is demonstrated that the telluric absorption band at 12.65 uw is a member of the 
family of carbon dioxide bands which determine the infrared limit of atmospheric 


transmission. 

Observations of the telluric spectrum made at the University of 
Michigan several years ago disclosed a moderately intense absorp- 
tion band in the neighborhood of 12.5 wu." Recent measures at the 
Lowell Observatory confirm the existence of the band and provide 
the position 12.65 u for its center.” 

The rather abrupt cessation of atmospheric transmission indi- 
cated by the Michigan observations at 13.5 4 has also been con- 
firmed by the Lowell measures. A detailed analysis of this cutoff 
now in progress at this observatory has revealed the important role 
played by the upper-stage transitions: 


42 


2v2\ . 4 32 , 342 _ k 
(v2); ( Y; ), ’ ( Y; ) Get 5 ’ Bes e: 2v2 + V; . 
2V; 


0 


NO 


As a consequence of the near commensurability of the frequencies 
v, and 2v, these upper-stage bands are widely dispersed on both sides 
of the fundamental frequency v,. This suggests the possibility that 
the band gq, is a high-frequency member of this group. An examina- 
tion of the energy relations of the carbon dioxide molecule or of the 
appropriate portion of the energy-level diagram (Fig. 1) taken from 


« A. Adel, V. M. Slipher, and E. F. Barker, Phys. Rev., 47, 580, 1935. 
2 A. Adel and C. O. Lampland, Ap. J., 87, 198, 1938. 


* A small and irregular trickle of energy appears to extend for a short distance be- 
yond 13.5 u. This behavior is a function of the carbon dioxide in our atmosphere. A 
detailed account will appear in a future paper. 
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the pioneer work of Martin and Barker proves that this is indeed the 


case.° 
It is possible that the weak telluric band g,, which has its center at 
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approximately 11.6, is similarly due to the highest frequency 


transition of the set 


f 4M 
5V2 
( ; ) — | 2v. + »; 
aT Pies , 
2V; / 0 


This point is being investigated. 
ARTHUR ADEL 
LOWELL OBSERVATORY 
FLAGSTAFF, ARIZONA 
April 1, 1938 


A NEW SPECTRUM VARIABLE: 5 LACERTAE 


ABSTRACT 

Sharp lines of H and K of Catt, superposed on the broad H and K lines of the 
cK5 component of the composite star, 5 Lacertae (K5Aos), are found to vary in in- 
tensity and radial velocity. The velocity of the late star, regarded as constant for many 
years, is also found to be variable. 

In the course of an examination of stars with composite spectra, 
the writer found that on a plate of 5 Lacertae taken in November, 
1934, the H and K lines of Ca 11 were of almost interstellar sharpness 
and were stronger than neighboring iron lines. They were super- 
posed on the very wide and shallow H and K lines which are char- 
acteristic of many composite spectra. Yet, on a plate taken in Octo- 


3 Phys. Rev., 41, 291, 1932. 
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ber, 1936, these sharp lines had completely disappeared, along with 
a few fainter lines in the K region, though the rest of the spectrum 
remained sensibly unchanged. By November, 1937, the sharp K line 
had definitely reappeared, together with a trace of the sharp H line 
and many faint lines in the wings of the broad K line. 

The star 5 Lacertae (a 22525™4; 6 47°11’—1900; m = 4.61) is 
HD 213310-1. It has been classed as KoAo at Harvard, Mo at 
Mount Wilson, and K2 at Victoria. Dr. Morgan has placed it at 
Ks5 on the Yerkes system. Its trigonometric and spectroscopic paral- 
laxes agree in assigning to it an absolute magnitude of — 2 or —3. 

TABLE 1 
RADIAL VELOCITY IN KM/SEC 














Date Ks Star Sharp K Line 
, Prob. 
Vel. é 
Error 
Nov. 5, 1934 +21.2 2 ee ee K, +21.4 km/sec 
WEES Boe BORO aieys ao arscnecats ts.2 1.8 No K or H 
Nov. 2 27 ae 16 I K, —20.0\,- 
~iiianieiy e 3 3 ; K weak; H, trace 
PHO 9S TOSI SS csce gece asses 12.2 1:2 K, —10.1 
WANED, BOSS «3 5ccs.6 are eo sas'4 17.0 0.6 
May 18 Se ei, Gio 5 0.5 Ee : 
May ag a = . : : ‘Spectrum does not cover region 
pA 27; .) Sir ra oe 
MMOL EORG <0 vee vs betes + 5.6 +5079 








In color it is redder than-any normal M star. The Henry Draper 
Catalogue states: ‘“From H to the end of shorter wave length, only 
the spectrum of the Ao star is visible.’ The star is listed as constant 
in light. 

The star is obviously binary, though no variable radial velocity 
has been reported previously. The radial velocity listed in the Lick 
catalogue of radial velocities’ is —11.6 + .o4 km/sec. The present 
additional measures indicate definite variability, presumably of long 
range. The radial velocities, reduced to the sun, for eight plates 
are listed in Table 1. The plate of January 19, 1938, was taken at 
the Yerkes Observatory and was made available through the 
courtesy of Dr. Struve. The last three plates were taken, and meas- 
ured, through the kindness of Dr. Petrie at Victoria, at the writer’s 


t Pub. Lick Obs., 18. 
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request. The remaining plates were taken with the Yerkes auto- 
collimating, one-prism. spectrograph attached to the Perkins 69-inch 
telescope. Rapid emulsions were necessary in order to get the inter- 
esting violet region. The Victoria, Perkins, and Yerkes plates had 
dispersions of 15, 30, and 55 A/mm, respectively. With the ex- 
ception of the Yerkes plate, there appears to be a systematic change 
of velocity with time, but this is not certain. The velocity varia- 
tion, however, can be considered well established. Through the cour- 
tesy of Dr. Moore, the original Lick plates of 5 Lacertae were re- 
measured and the previous value of —11.6 was confirmed. 

The sharp K line does not show the same radial velocity as the 
lines of the late-type star, but apparently it has a velocity variation 
of its own and, thus, is not interstellar but may belong to the early 
companion. It is most interesting that the sharp H and K lines were 
strongest when their radial velocities were the same as that shown 
by the late star, and were much fainter when their velocities differed 
from those of the late star. At the latter epoch the broad H and K 
lines of the K5 star also seemed stronger. On October 25, 1936, no 
sharp H and K lines appeared and the radial velocity of the late- 
type star was nearly the same as on November 5, 1934, at which time 
the sharp H and K lines and the late-type lines had the same radial 
velocity; hence, the secondary may be partially or totally eclipsed. 
However, the constancy of light seems to preclude this explanation. 
Artificial composite spectra made directly at the Perkins telescope 
have shown that in combining a Ko and an Aos star, the resulting 
composite K line always shows a sharp core, and it is difficult to 
explain the disappearance of the early K line without an eclipse, 
thus leading to an apparent contradiction. More observations in the 
near future are necessary before any definite mechanism can be as- 
sumed. 

The early star is classified correctly as Aos, since the narrow H and 
K lines could not belong to a later type, while no B-type lines, espe- 
cially He 1 4026, are found. The late component shows no 77 O 
bands in the green, and hence the Victoria K2, or the Yerkes Ks, 
classification can be adopted. The star is clearly a supergiant ac- 
cording to both its color and its parallax. 

The writer hopes that the star will be placed on a photometric 
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program at some observatory to check for the small light variation 
which can be expected if the spectrum changes are the result of an 
eclipse. 
J. A. HYNEK 
PERKINS OBSERVATORY 
DELAWARE, OHIO 
June 15, 1938 





ERUPTIVE PROMINENCES AND IONOSPHERIC 
DISTURBANCES 

It has recently been shown' that every solar eruption causes a 
disturbance in the ionosphere, and that the intensity of the disturb- 
ance does not depend on the distance of the eruption from the center 
of the sun’s disk.? The case of eruptions occurring within a few de- 
grees of the limb is presented in the accompanying table, in which 
important eruptive prominences observed with a spectrohelioscope 
at the Commonwealth Solar Observatory, Canberra, are given, to- 
gether with ionospheric effects, as determined by Mr. A. J. Higgs 
with the Observatory P’f equipment. The estimated central meridi- 
an distance is given for the base of the prominence. 














Date E.AS.T. (ie. | = seeeeaeee 
Kilometers Disturbances 
‘ 
1936Sept. 2.. ght5™—10'30™ 82 225X103} F, disturbance 9"15"™—9"30™ 
1937 July 3...| 13 28 —13 49 go 75 No effect 
20: ..| 26:02 16 38 >90 350 No effect 
39'<..-4/ SOLERO II 22 >9go 120 No effect 
Aug. 27...] 14 08 —14 58 85 60 Fade-out 13"10™—13"50"™ 
Sept. 15. iI —- 309 >go 225 No effect 
iS 12 22 —I2 44 >9go 75 No effect 
Get. 22:..1 12 36 iz 2 85 210 Fade-out 12'30™—13"10™ 
25...| 9 10 —II 22 go 60 Fade-out 8530™—g"35™ 
Dec. 9...] 15 39 —16 34 >9o0 235 No effect 
1938 Feb. 24...] 12 49 I2 53 >go 150 No effect 
Mar. 16...] 12 12 13 04 >go 110 No effect 


From the accompanying table it is seen that no solar eruption the 
base of which is behind the limb causes an ionospheric disturbance. 
It would seem then that La or other ultraviolet radiation which can 


* Nature, 140, 603, 1937. 2M.N., 97, 605, 1937. 
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affect the ionosphere has its origin at the seat of the eruption, and 
that it is not emitted in appreciable quantities from the bright ma- 
terial comprising the eruptive prominence itself. 

R. G. GIOVANELLI 


COMMONWEALTH SOLAR OBSERVATORY 
CANBERRA, AUSTRALIA 
March 24, 1938 





VARIABLE HYDROGEN EMISSION IN THE 
SPECTRUM OF y URSAE MAJORIS 

y Ursae Majoris (Aon) shows not only Ha emission but a 
rather peculiar spectrum for a star of its class. A series of plates of 
the spectrum in the visual region recently taken with the Yerkes 
auto-collimating spectrograph attached to the 69-inch reflector of 
the Perkins Observatory reveal a weak emission line in the core of 
the broad Ha absorption line. The star is not listed in Merrill’s 
“Catalogue of Be and Ae Stars.””" 

Our first plate, taken April 19, 1938, on Ha emulsion, showed 
definite but weak emission centered in the core of Ha. The absorp- 
tion profile was unsymmetrical, being steeper on the violet side, so 
that the emission component appeared to lie on the violet side of 
the center of the line as a whole. The absorption line has extensive 
wings which cover an interval of about 60 A. The core appears to 
be about 1o A wide. On April 27 (III F plate) the emission com- 
ponent was much less conspicuous, and it appeared to be shifted 
to the violet edge of the core. Our spectrogram of May 3 showed 
essentially the same line structure. On the following night the emis- 
sion component was even more diffuse. In fact, there was some 
doubt whether it actually existed. A microphotometer tracing sug- 
gested its presence by a slight hump on the violet edge of the core. 

An attempt was made to secure another photograph on III F 
emulsion on May 23, but clouds cut it short. The spectrum was so 


tAp. J., 78, 87, 1933; Mt. W. Conir., No. 471. 
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thin that little could be seen, and a microphotometer tracing showed 
nothing definite other than slightly greater strength in the violet 
half of the profile. On June 5, however, the emission core was again 
definitely shown. It appeared as strong as it had on our first plate 
and considerably stronger than it had on any of the later plates. 
This time it again appeared centered on the absorption core, but 
slightly to the red of the center of the line as a whole. Our plate of 
June 9 showed no emission component at Ha. 

The rest of the spectrum was investigated. A preliminary glance 
showed nothing but the two hydrogen lines, but a careful examina- 
tion revealed a number of very wide and extremely shallow absorp- 
tion lines. They were from 1o to 15 A wide and were so shallow as 
to be almost invisible through lack of contrast. Sixteen lines were 
definitely located, and several other fainter ones are probably 
present. Although only the roughest measurement was possible, 
the lines were identified with those of ionized metallic atoms. The 
consistency of the correspondence between these shallow absorption 
lines and the strongest emission lines in the spectrum of y Cas- 
siopeiae lends plausibility to the identifications. The accompanying 
list was found. The last line was much sharper and darker than the 


ee ee Feit PMR le ys ire tes Fei 
accents Fe ir Md 2 rach nasa nla ee Nat 
BE rrcan uy eae SY 11 5896 Nal 
Na iss 2 bes Feu SN rip oa Osan s 4 Si II 
meee Ti eee $7 II 
es eee Feu ee Si 11 
Lee catsorn Ti ee coreces oben Feu 
A 288 eA Fei | eee Fei 


others. The Na D lines were extremely diffuse. Some of these lines 
were recorded by Miss Waterman.’ 

On the plates of June 5 and 9 a weak emission component ap- 
peared in the core of HB. Owing to curvature of field in the grating 
spectrograph the focus is never sharp at HB. Consequently I do 


2 Lick Obs. Bull., 8, 1, 1913. 
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not attach much weight to the appearance of that line on any plate. 
I merely mention this observation for the interest of others who may 
care to follow the star with a prism spectrograph. 
ERNEST CHERRINGTON, JR. 
PERKINS OBSERVATORY 
DELAWARE, OHIO 


June 9, 1938 


REMARKS ON THE PAPER “THE SODIUM CONTENT 
OF THE HEAD OF THE GREAT DAYLIGHT 
COMET SKJELLERUP 1927 K” 

In Section IV of this paper we have estimated the magnitude of 
the cometary head in the light of the D lines, combining the ap- 
proximate dimensions of the comet’s atmosphere with the density 
of the normal sodium as determined in Section III. In doing so we 
assumed the monochromatic magnitude to be proportional to the 
mean density. But according to the foregoing considerations in our 
paper, especially Figure 4, page 356, this assumption is only ful- 
filled up to C ~ 0.2 or Nf ~ 3.3, where N is the number of normal 
atoms per cubic centimeter and f the line strength. Putting f = 1 
and N = 3.3 for the sum of both D lines, and using the values 
given on pages 357-358 of our paper, one gets for the ratio of the 
sun’s apparent brightness and that of comet ‘ kjellerup’s head F = 
3.5 *% 20". 

Our Figure 4 quoted above shows further that the intensity of the 
D lines increases 4.5-fold if NV increases from 3.3 to 50—that is, to 
its lower limit calculated in Section III. With V = 50 one gets 

f= 2 Ye 7.8 X 10", 
4.5 
and putting this equal to 2.5”, « = 27.4; therefore the temporary 
monochromatic magnitude of the comet’s head was 


M = 27.4 — 26.7 = +0.7, 
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a value which is compatible with the subjective estimate of Mr. 
Slipher mentioned in Section I. 

The upper limit of the sodium density found in Section III, 
N = 2500 (C = 150), would give approximately the same intensity 
of the D lines (see Fig. 4)—so that under the peculiar conditions of 
the comet the estimated monochromatic magnitude of the comet’s 
head does not allow us to decide between the upper and the lower 
limits of the amount of sodium calculated in Section ITI. 

ARTHUR ADEL 

V. M. SLIPHER 

R. LADENBURG 
LOWELL OBSERVATORY, FLAGSTAFF, ARIZONA 


PALMER PHysICAL LABORATORY, PRINCETON UNIVERSITY 
June 1938 





ERRATUM 
Dr. T. E. Sterne calls attention to a misprint on page 149, section 5, line 10, 
of his article, “The Secondary Variation of Delta Scuti,” in the March, 1938, 
issue of the Astrophysical Journal (Vol. 87). The word “possible” should read 
“impossible.” 











